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FLAT PLATE CASCADES AT SUPERSONIC SPEED* 
By Rashad M. El Badrawy 


INTRODUCTION 


The cascade problem in the subsonic range can be analyzed under 
certain assumptions either by mapping or substitution of the blades by 
singularities - sources, sinks and bound vortices - where the separation 
of flow from the blades can cause various departures from the obtained 
results . 

Raising the flow velocity to a given value is accompanied by sonic 
velocity within the cascade, which usually renders the solution of the 
problem even more difficult. The same complication exists on the cascade 
in flow at supersonic speed, in which the velocity is retarded to sub- 
sonic by shocks. 

But when the cascade operates entirely in the supersonic range, the 
conditions become clearer. All disturbances act downstream only from the 
sources of disturbance, so that the pressures and velocities at the sur- 
face of a sufficiently thin airfoil in the stream can be readily determined. 

The present report deals exclusively with problems of cascade flow 
in the supersonic range. As is known the flat infinitely thin plate is 
the best airfoil with respect to wave resistance in supersonic flows; 
hence it is logical to start with the cascade of flat plates. The last 
chapter deals with the case of finite thickness. 

Lift and wave resistance of an isolated plate are computed first 
since the cascade problem can often be reduced to this special case. 

The well-known theories of two-dimensional supersonic flow are applied - 
that is, the laws of oblique compression shocks and the expansion around 
a corner. 

The air forces are then calculated again and compared with the pre- 
viously obtained exact values by means of Ackeret 1 s formulas of linearized 
theory . 


*"Ebene Plattengitter bei Uberschallgeschwindigkeit . " Mitteilungen 
aus dem Institut fur Aerodynamik an der E.T.H. , no. 19, 1952- 
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The cascade problem was to be solved in such a way as to be free 
from the inevitable inaccuracies of the graphical method. For this reason 
the cases of overtaking, crossing and reflection of compression shocks and 
expansion waves frequently occurring on supersonic cascade flows, which 
usually are solved by graphical method, are analyzed in chapter II. 

In chapter III the cascade problem is discussed and its solution 
described in the light of the results obtained in chapter II. A numerical 
example is also given. The same chapter gives further a definition of 
the efficiency of the simple supersonic cascade and an evaluation for 
several angles of stagger and attack. 

The small angles of attack involved justified the use of a linearized 
cascade theory. 1 This is done in chapter IV. The numerical example of 
chapter III is thus linearized and the results compared with those of the 
exact solution. The supersonic cascade flow at various angles of attack 
was recorded by schlieren photographs of the flow between two parallel 
plates, in the high-speed wind tunnel of the Institute (chapter V). 

Chapter VI deals with the specific case of unsteady flow through 
the cascade, caused by abrupt angle-of-attack changes. 


^According to Ackert's linearized theory, the lift and drag of a 
double-wedge profile of thickness d and chord l at angle \|< in super- 
sonic flow M is, in the presence of friction (cp) 


the wave resistance should be equal to the sum of friction drag and thick- 
ness effect. In that event 



For the best drag-lift ratio e = — , put = 

c a 


0. This means that 


c 



\Jvfi - 1 

4 



Assuming possible values for d / l and cp results in comparatively 
small optimum angles i Q pt‘ 
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In chapter VTI the effect of friction and thickness in a special 
case on the cascade efficiency is analyzed. Since there might be a 
possible application of the supersonic cascade to the supersonic pro- 
peller, a simple evaluation of the efficiency of such a propeller is 
made. A parallel steady two-dimensional flow is - with exception of 
chapter VI - postulated. 

The conventional notation is used unless specifically stated other 
wise in the text. 


If 
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CHAPTER I. THE FLAT PLATE 
1. General Considerations - Stipulations 


The general equation of continuity of any compressible flow is 


Sp t d(pu) ( d(pv) t a(pw) _ Q 

St Sx Sy Sz 

The rate of propagation of a small disturbance, that is, the sonic 
velocity, is, as is known 


a 4 - 


Sp p 

r— • = K — 

Sp P 


( 2 ) 


C P 

where k = — . 

c v 

In flows, in which a flow potential cp exists, the continuity equa- 
tion can be written as 


C \ 2 2 2 

+ IfP + IjP . It = 0 

at a x a x a y a y a z a z y dx 2 ^2 dz 2 

The momentum theorem gives the following relations 


1 ap 

^9 + * 2 . A + 

M A + 

acp a 2 cp 

p a x " 

ax at ax ^ x 2 

a y ax a y 

a z ax a z 

1 ap _ 

a 2 cp acp a 2 cp 


acp a 2 cp 

p a y 

Sy St Sx Sy Sx 

* 3y 2 

a z a y a z 

1 ap 

2 2 
S cp J Sep S cp 


+ §9. A 

p a z 

Sz St Sx Sz Sx 

a y a z a y 

a z 2 


(3) 


( 9 ) 
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In two-dimensional flow the potential must therefore satisfy the 
equation 



/ \ 2 

p 

/ \ 2 

1 - -i/Vl 


1 - 

a 2 Wj 

ay 2 

[ a 2 W 


2 Sep Sp S cp 2 Sep S cp 

a 2 Sx Sy dx dy a 2 Sx Sx St 


2 2 

2 Sep S ep 1 S cp __ Q 

a 2 Sy Sy St a 2 St 2 


(5) 


The velocity of sound is then 

a 2 = a ^ 2 - (k - 

where aQ = velocity of sound 
For the steady case, the 



, >2 

, \ 2 

\ 

1 

m + m 

+ & 

2 

k. 

V*J 

\8y J 

St 

> 


in state of rest, 
equation is reduced to 


(6) 


ax 2 



0 

. \2 

Ns)J 

. djP 

1 - l/M 

8y 2 

a8 WJ 


^2 

a cp 


2 a? acp 

a 2 ax ay ax ay 


= 0 


( 7 ) 


and the flow is completely identified, if the function cp(x,y), which 
is to satisfy the boundary conditions, is determined. 


This equation is either elliptic, parabolic, or hyperbolic, depending 


upon 


(l - M 2 ) £ 0 


where 



( 8 ) 


is the local Mach number. 
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The use of this equation is difficult if its type in the particular 
range, as in the transonic range, is changed. However, the flows analyzed 
here, are of identical character everywhere, that is, the flow is of the 
hyperbolic type. 

One of the known solutions is the expansion around a corner, developed 
by Prandtl and Meyer (ref. 1) . 


The two-dimensional flow past the wall AE at a Mach number Mq 

(fig. 1) is deflected by a convex bend at E through an angle ©, through 
which an expansion is initiated. The disturbance proceeding from E 
spreads out solely in the range lying downstream of the Mach line EBq, 

where 


In it M 2 is the Mach number of the flow after the expansion. 

The streamlines in range BqEB2 are curved similarly and run 
parallel to the wall ED downstream of this range. 

It can be proved that the Mach lines in this flow are the character- 
istics of the differential equation which define the potential. 


2. Conditions at Expansion Around a Corner 



(9) 


and stops at the Mach line EBp, where 



When the expansion proceeds from a Mach number Mq = 1, 
the following relations can be proved (ref. 2) : 



tan p -2 = A cot Aui 



( 10 ) 
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P2 

PO 


' K + 1 ^ 

2 cos P'Aoii 


K 

K-l 


(Pq = stagnation pressure) 


Obviously 


M 2 =' 


(k + 1) - 2 cos^Ao) 
(k - l)cosP>va> 


0 


a> + b 2 


2L 

2 


As function of M 2 (ref. 3) 


( 11 ) 


( 12 ) 


(13) 


0 = 


cos -1 — + — A cos~l / 1 - 
Mo 2 


K + 1 


1 + 


( 14 ) 


This equation gives a maximum angle of expansion ©max* which cor- 
responds to a Mach number M 2 = ~ after expansion (k = 1.400) 

©max = 130 • 45° 


If the Mach number before the expansion Mq is assumed other than 1, 
the maximum angle of expansion becomes obviously 


©max Mq " ©max ” v 

where v is the angle of expansion from M = 1 to Mq. The values for 
various Mach numbers of the inflow (Mq) are 


Mq 

1.00 

1.50 

2.00 

2.50 

5 

8 

10 

00 

®max Mq° 

130.45 

118.55 

104.07 

91-32 

53.35 

34.53 

28.14 

0 
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3 . Conditions of Oblique Compression Shock 

The discontinuities that may appear in supersonic flows and across 
which velocity, pressure, density, temperature, and entropy undergo a 
discontinuity, while the total energy, thermic and mechanical, remains 
constant, were predicted by Riemann (i860) and Rankine and Hugoniot (1887) 
as normal compression shocks. 

In oblique shocks (Prandtl-Meyer) only the velocity component normal 
to the shock front is modified. 

In figure 2 the supersonic flow past the wall AE is deflected at E 
by an angle 5. A compression shock is produced and the shock front ES 
is inclined at an angle 7 - the shock angle - toward the air flow 
direction. 

With subscript 1 denoting the state before the shock and subscript 2 
that after the shock it can be proved that (refs. 3 and. 4) 


-2k/ m 2 Bin 2 7 . 

P 2 _ k + 1 \ 1 2k J 


( 15 ) 


1 
1 K 


K - l/Pl, , 


4k /* - 1 

sin ^7 


P 0 is + 1\P 0 J ( K + l ) 2 \R + 1 


1 + 


bn. . 

sm^7 


(k - l ) 2 


cot 5 = 


2 ( 

( 1 + 

IS + ll 

,M2_ 2 sin 2 7 

: + 1 

CM 

H 

X 

2 M]_ 2 sin 2 7 - 1 

V = 

cos 7 

u l" 

cos (7 - 5) 


1 tan 7 


K — 1 

'M K 


(16) 

( 17 ) 




(18) 


( 19 ) 
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tan(y - 8) p l 1 


tan 6 

^ 2 ' K + 1 

1 K - 1 

-1 


2 

1 

Mg^sin^y - 5 ) 


( 20 , 21 ) 


A direct relation between the Mach numbers before and after the 
shock can be established 


^2 cos y | j 

/ p 2 '/ 

1 Pi 

cos (7 - 5) u 


/Pi 


( 22 ) 


The relation for the change of the static pressure by the shock is 
the same as for the normal shock when it is applied to the velocity com- 
ponent perpendicular to the shock front. Consequently 



rc + 1 p . p 
siry 


V 3 


1 + 


— M-^sin^y 


f — — Mn 2 sin 2 7 

Ik + 1 1 


7 



(23) 


From these equations it follows that the shock angle 7 is greater 
than the Mach angle, that is, the speed of propagation of a finite dis- 
turbance is greater than the sonic velocity. When the angle of deflec- 
tion 5 approaches zero, 7 = p and the shock changes to a Mach wave. 

Also of interest is the shock angle at which the Mach number after 
the minimum shock becomes equal to unity. Denoting this angle by 7 S 

it can be proved (the weak stable compression shock is always allowed 
for) (ref. 3 , P- 47 ) that 



( 24 ) 
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This equation is used to determine the maximum shock angle which 
corresponds to a Mach number before the shock Mq = » and a Mach num- 
ber M2' = 1 after the shock. The result is 

sin 2 7 s = K - + 1 (25) 
2k 


hence 


7 S = 67.8° 


at k = 1.400 (air). 

By equation (18) the corresponding deflection angle 5 S is 


6 S = 45.58° 

O 

Table 1 and figure 5 represent the values of 7 S and & s at 
various Mach numbers Mq. 

4. Lift and Drag of an Infinitely Thin Plate (Exact Solution) 

An infinitely thin plate ab in parallel flow at supersonic veloc- 
ity Uj_ is placed at the angle \|r. It is assumed that the width of the 

plate transverse to the flow direction is «>, so that the problem is two 
dimensional. 

The streamlines above oa (fig. 4) experience a deflection which 
is associated with an expansion. So the state at the upper side of the 
plate can be defined by equations (10) to (l4). But below the plate a 
compression shock ad occurs. The state of the flow on the lower side 
of the plate is accordingly determined from the formulas (l6) to (21). 

The force on the plate per unit area is 

K = (p 2 ' - P 2 ) (26) 


The weak stable shock is always taken into account. 
ZAMM, 19^8 and Thomas , Proc. N.A. Sc., Nov. 1948. 


See Richter, 
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where P 2 ' and. P 2 represent the pressure on the lower and upper side 
of the plate. 

Obviously, the lift A and the drag W per unit width are 


A = K cos i|/L 

(27) 

W = K sin tyL 

( 28 ) 


To compute a lift coefficient, a reference dynamic pressure of the 
inflow 


1 2 
*L * 2 PlUl 

or 


<11 = 


K 

2 


Pl M l 2 


(29) 


is utilized. 

As function of the Mach number the ratio of dynamic to air stream 

pressure is 


pi 



that of dynamic to stagnation pressure is 


p 0 


K 1 M 2 
M-, = 

2 Po 1 



K 



( 30 ) 


The results are represented in table 1 and figure 5* 
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Lift and drag coefficients are herewith 


A p 2 ' p 2 
c a = — — = cos \|r 


c w ~ 


W / p 2 ' " P2\ . 


*L l 


sin >(/ 


or, if all pressures are referred to stagnation pressure Pqj 


c a - 


P2 P 2 
& ' *>. 


«i 


COS \Jf 


c w = 


p 2 p 2 

jo~ " 51 


* 1 


sin i|r 


(3D 


(32) 


The drag/lift ratio is 


e = — = tan 


(33) 


Table 2 gives the values of c a , c w , and € up to Mp = 10 as 
computed by the formulas ( 32 ) and (33)* 

In the calculation of the Mach numbers up to Mp = the tables 
by Keenan and Kaye (ref. 6) as well as those by Ferri (ref. 3) were used 
to define P 2 /P 0 an< ^ P2/Pl ( K = 1-^00). 

For higher Mach numbers, the formulas of sections 2 and 3 were 
employed. At each Mach number, the angle of attack was varied up to 

* s (m 2 ' = 1 ). 

Figures 6 and J show the variation of c v and c a over the angle 
of attack figure 8 shows the polar s c a plotted against c w . 

The boundary curves show the maximum lift and drag coefficients that 
can be expected without getting in the transonic range. 
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Other values for the boundary curve are given in table 5* Since the 
pressure distribution on the upper and lower side is constant, the result- 
ant force is applied at plate center and is normal to the plate. There is 
no suction force as in subsonic flow. 


5. Lift and Drag at High Mach Numbers 

At high Mach numbers the angle of attack of the plate can exceed 
the maximum expansion angle ®rriax (section 2) corresponding to the Mach 

number of the airstream Gfoiax = M]_ S' 6.4). Hence, when assuming 

continuous flow, an empty wedge-shaped zone between plate and flow appears. 
This zone is largest at constant angle of attack when M^ = In that 

event, no deflection of flow is possible. 

Owing to this vacuum space, the pressure at the upper side is zero. 

The resultant force K is obtained then from the pressure on the lower 
side, behind the compression shock. Hence, per unit area 


K = p 2 ' 

or, when referred to the dynamic pressure of the airstream, 


(34) 


^ I l 

K_ _ P2_ _ 2 £2_ 

*1 ' q l = kM;l 2 P 1 


(35) 


Introducing P2VP1 from equation (15) gives 


'll 


b . 2 

sm^7 

k + 1 


k - 1 2 


K + 1 



where the term containing 


1/M1 2 


can be disregarded without great error. 


K 4.0 

— = sir7 

<ll « + 1 


( 36 ) 
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So the lift and drag coefficient are 


c 


a 


°w 


k , p 

— — — sin ^7 cos 
k + 1 


k p 

— - — siry sin \Jf 
k + 1 


(37) 


Both formulas are dependent on 7 and only. Between these there 
exist the relation given in equation (l 8 ), which can be written as 
follows (6 = \(r) : 


cot 


f K + 1 
2 


1 


• 2 1 
sin^7 - — - 


Mi 4 



7 


where the term 



can be disregarded again. 


Then 


cot 


I K. + 1 
t 2 sin^y 


l\ tan 7 


(38) 


The values and curves designated with M = «> in table 2 and figures 6, 

7 and 8 were defined by equations (37) and (38). For comparison the lift 
and drag was also computed by Newton's formula (the normal component) 


Ca = 2 sin^ty cos 
c w = 2 sin3\|/ 


I 


( 39 ) 


The corresponding values and curves carry the subscript N. 
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6. Calculation of Lift and Drag by Linearized Theory 
According to Ackert ' s linearized theory, the members of higher order 
in — can be disregarded without great error in the potential equation 

dy 



1 - -M 2 


1 - 

a 2\bxJ 

By 2 

a2\dy/ 


2 

2 &p dcp 5 cp 
a 2 dx dy dx dy 


for slender bodies at small angles of attack, because the interference 
flows are small compared to that of the air stream. 

The equation reads accordingly 




+ 0 

dy 2 


Inserting 



and observing that M is greater than unity, the equation reads 


. d . is. , 0 (40) 

ox 2 dy 2 


The general solution of this equation is 






( 41 ) 


It indicates, as stated in section 2, that the lines of constant 
potential are the Mach lines of flow, and their slope has the Mach angle \i. 


l6 
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This solution shows further that the flow velocities 


Sx 


u = — and 


dcp 

V = — 1 - 

dy 


satisfy the condition 


u = 


\i 


( 42 ) 


NF - 1 


At the surface of a body in the stream 


u = 

U dx 


is applicable. 

The pressure variation by the momentum equation reads 


— = - U AU = - Uu (A 3 ) 

P 


where U is flow velocity and u is interference flow in stream direc- 
tion. 

Accordingly 





M 


dx 


1 


( 44 ) 
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The pressure difference between both sides of a flat plate is 


£p 4 dy 

q |/ m 2 - 1 ^ 


(45) 


The lift and drag coefficients at the angles in question are 






M 2 - 1 


C W - 


kr 


/ 


M* - 1 


The drag/lift ratio according to this theory is 


(46) 


e = 



(4?) 


Instead of the expansion wave and the compression shock at the 
leading edge, it has simple Mach lines as interference lines (fig. 9), 
in contrast to the exact theory. 

The values given in table 4 and plotted in figures 10 and 11 were 
computed by these formulas. The calculations were carried out at each 
Mach number up to angle of attack \|j s - from the exact theory. The cor- 
responding c a and c w values lie on the curve G f . 

At sonic velocity on the lower side of the plate P2VP0 = O.5283. 
This value, introduced in the following directly obtainable relation 




1 


( 48 ) 
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and. 


P 2 7p 0 - Pi/pp _ 
^l/Po 



1 


gives the angle of attack (i s j) corresponding to Mp' = 1, which in gen- 
eral is greater than ty s * 


7. Comparison of the Results of the Linearized Theory With 
Those of the Exact Method 

In table 5> "the difference is (cq - c^) , where cq is the coef- 
ficient of the exact method and c L is that of the linearized theory 
at Mach numbers Mp = 1.40 and Mp = 5 >00. 

It follows that the linearized theory is a very good approximation 
for small angles (up to about 10°) • For greater angles the values of c a 

and c w are too small. 

In figures 6 and 7 , the c a and c w curves by linearized theory 
marked A' and A are included for comparison. 
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CHAPTER II. INTERSECTION, OVERTAKING AND REFLECTION 


OF COMPRESSION SHOCKS AND EXPANSION WAVES 


1. Introduction 


Overtaking of expansion waves and compression shocks in supersonic 
flows occurs when the marginal streamlines - or boundary walls - change 
their direction twice in the opposite sense (fig. 12(a)). 

If expansion waves or compression shocks strike a fixed wall and 
their slope toward the wall does not exceed a given angle, they are 
reflected as expansion waves or compression shocks (fig. 12(c)). Crossings 
occur in flows through channels and free jets (fig. 12(b)). All these 
events can occur in cascade flows (fig. 12(d)). 


(a) Suppose that a small expansion occurs at B in the supersonic 
flow Mi, pi, ap, past the wall AB (fig. 13)- The angle of 

expansion is ZB. If ZB is sufficiently small, differential considera- 
tions are permissible. 

Bernoulli f s equation gives 


where U is the magnitude of the velocity and AU its variation; Ap is 
the pressure variation. 

Since the vectorial velocity variation is normal to the interference 
line, the variation of U is 


2. Small Variations 


£Pl = - PjU-l 


(49) 


a| U | = Uj tan A0 = 


!!]_ A0 


(50) 
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hence 


z® = 


P^ 2 ^© 

\j^l 2 - 1 


But, as the dynamic pressure q is given by 


(51) 


q. = 



the pressure variation can be written as 


■ - 1 - - - - (52) 

- 1 

M follows at (ref. 3 > P* 26 ) 

/ k - 1 

_ -1 + M 

\Ju? - 1 \ 2 

The Mach line BE^ forms with flow direction AB the Mach angle 
the Mach line BE2 at the end of the expansion the Mach line P2« Now 

it may be assumed that this small expansion takes place on the inter- 
ference line BE ' , whereby 



= 


The variation of the Mach number 


Mn 


£ A'BE 1 


p x + p 2 " A0 
2 


( 54 ) 


(b) A simple differentiation gives the change of the shock angle 7 
as well as the pressure change P2* after the shock, due to a small vari- 
ation of the angle of deflection 5 , for the compression shock AB 

(fig. 14 ). 
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Between 7 and 6 the relation (eq. 18) 

“i 2 \ 


cot 5 = tan 7 


- 1 


- 1 


/ 


exists, and therefore 


where 


= - sin^5 


and 


1 

(a \ 

sec 7 

r 1 ) 


MjL^sin^7 — 

B 2 


A7 = c A 7 (55) 


. K + 1 2 

A = Mq^ 


B = ( Mq^sin 2 7 - lj 


C = - sin 2 & 


sec^yl 

'7 . A 

! 

V B ) 


2 • 2 2A 
L sin r — 
B 2 


The pressure p,-,' after the shock is (eq. 15) 


, , 2k ,,2-2 k + 1 

p ’ = Pl Mq sin^7 - 

1 K + 1 K — 1 


(56) 


and the result for a small variation of the shock intensity is 

, p 2k 

Ap 2 = PqMq c sin 2 y Ay 


k + 1 


( 57 ) 
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3- Overtaking of Compression Shock and Expansion Wave 


(a) The supersonic flow Mp, Pp, p ± past the wall AE (fig. 15(a)) 

undergoes a directional change 5 at E. The compression shock EF and 
the flow direction form the shock angle y 1 in zone (l). At C an expan- 
sion takes place about an angle 0, and the expansion wave FCG overtakes 
the compression shock at F. To simplify the calculation, the continuous 
expression is replaced by a given number of expansion waves of finite 
intensity, whereby a successive expansion through these waves is assured. 

If n is the number of waves, the expansion due to a wave is 0/n = A0. 

The number n must be so chosen that A© is sufficiently small. 

Now consider the intersecting of wave CF^A©) and compression 
shock EFp(b), figure 15 (b). 

From Fp the compression shock advances with weaker intensity in 
'ii^’sction Gp, that is, it deflects the flow less — say by 6' • FpGp forms 

with the flow direction the angle y at (l). Indicating the various zones 
by (l); (2), (3), and (4), the streamline through Fp splits the zone (4) 

into the portions (4 C ) and (4u). The flow in (2) and (3) is fully known, 
because the angles & and A0 are known. 

To define the conditions in (4), the streamlines Sp, S 2 , and S 3 

are examined. The directional change of S 2 amounts to (5 - A0) . But 
along S 3 the flow experiences the directional change 6 ’. To maintain 

equilibrium in (4), the pressure as well as the velocity direction above 
and below the streamline Sp must be equal. 

In general, the pressure change from (l) toward (3) is not the same 
as from (l) to (4), so that a reflected expansion wave - possibly a small 
compression too - must appear between (3) and (4), say along a line FpHp. 

Supposing that this reflected wave is an expansion wave of inten- 
sity A0 . By intensity' of an expansion wave or a compression shock is 
meant the deflection, which the flow experiences in the process . 

The pressure in (4u) is, (according to section 2) 



1 


( 58 ) 
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The difference of the deflection angles amounts to ( 6 ' - 6 ) = AS. In 
general, this intensity decrease is small, because the compression shock 
is much stronger than the expansion wave. The corresponding change in 
shock angle 7 is Ay. 

The pressure in (4 0 ) follows from the change in shock intensity. 
Hence we can say that 


Pl+o = P2' - 4?2* 


(59) 


This is again the equation for small variations derived for shocks from 
equation (57)- Accordingly 


o 2 ft 

p 4 0 = p 2 ~ p l M l Sin 27 A7 


(60) 


Posting p^ = pi, , gives 


p 5 




1 - 


i 


m 3 2 - 1 


A© j = p 2 ' - Pi M i 2 —— — sin 2y Ay ( 6 l) 


K + 1 


For the velocity direction in zone (4) to be unequivocal, it must 


AB = A 6 + £&' (62) 

The relation between shock intensity variation and angle of shock 
(eq. ( 55 )) together with the two previous equations gives 


A©' = 



I A©FG - A _ f 


(63) 
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A ££> A3 - A3' 

A 7 = — = 

C C 


(64) 


with the constants 


F = 




M-^sin 2y 


(u^sln 2 ? - 


The condition for equality of static pressures is not identical with 
that for equality of velocity magnitude above and below the streamline S^. 

As the shock losses on either side of the intersection point F are 
unlike, the stagnation pressures in the wake above and below streamline Sq 

are different, hence there is a small vortex layer along this streamline. 

Figure l6 represents the graphical solution of the problem by means 
of the characteristics and the shock polars. The condition for equality 
through equality of velocity magnitude in the entire zone (4) is 
approximated . 

(b) The reflected wave is disregarded: 

In general, the angle of deflection £&' - intensity of the reflected 
wave - is very small (compare numerical example). Thus the pressure in 
(3) is not much unlike that in (4), so that this reflected wave F 1 H 1 can 

be discounted. 

In this event the flow directions in zones ( 3 ) and (4) are identical, 
or in other words 


A0 = (65) 

With equation (55) A5 can defined, and from it the new direction y' 
of the compression shock. The velocities in (3) and (4) have then obvi- 
ously the same direction but not the same magnitude by reason of the small 
vortex layer developing between (3) and (4). 
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Free stream: 


Numerical Example 


P1/P0 = 0.12780 

Before overtaking: 

shock intensity 
shock angle 

hence the state in zone ( 2 ) 

P2/P0 = °‘ 1 ^ 

intensity of the expansion wave 
state in zone ( 3 ): 

P3/P0 = 0*1695 

Determination of constants: 


M x = 2.000 

5=6° 

7 = 35*24° 

m 2 = 1.780 

A 0 = 1 ° 

M5 = 1.818 

c = 1.078 
F = 3*04 
G = 3-015 
Q = 3*175 


Inserted in equation ( 63 ) and (64) gives: intensity of the reflected wave 


By equation (62) 


Ay = O.89 0 
A0' = 0.06° 

A5 = A0 - A0' = 0-94° 


Therefore the shock intensity after overtaking is 

6 = 5-04° 

The new shock angle is 

7’ = 7 - A7 = 34-35° 

The reflected wave disregarded, leaves 

A 7 = C = C A0 = I.O 78 0 

shock intensity 5 1 = 5 ° 

angle of shock 7 ' = 34.16° 


It is readily apparent that the reflected wave is very small, hence 
scarcely affects the pressure in zone ( 3 )* 


2 6 
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4. Intersection of Compression Shock and Expansion Wave 


Figure 18(a) represents an expansion wave AG of intensity 0 
proceeding from the corner A. At F this wave crosses a compression 
shock of intensity 8 emanating from the corner B. 

As before, the continued expansion is replaced again by n expansion 
waves between which the flow is straight. The deflection by each wave is 


After crossing (fig. 18(b)) the compression shock has an intensity 8' 
and a shock angle y' . Now the expansion wave has the intensity 20 ' . 

The zones produced this way are numbered (l), (2), (3), and- (4). 

The streamline FpS splits the zone (4) into (4 0 ) and (4 U ). 

Looked for now is the shock intensity &', shock angle y l , and 
expansion angle 20' after crossing, and the state of flow in (4), when 
the state of flow in (l) 8, 7 and 20 are known. 

According to chapter I the state of flow in (2) and (3) can be 
determined directly. 

The pressure in (4^) follows through a small expansion 20' according 
to the laws in chapter II at 


where A0 is still unknown. 

The method of solution consists in first making an assumption for 
the shock intensity after intersecting, which is 



(66) 


= (8 - 20 ) 


(67) 


The corresponding angle of shock would be 7 i - 
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The state after this shock, indicated by b Q ' , can again be defined 
according to chapter II. The pressure in (4 0 ') is 


p 4o 


2k 


Mj; sm^. - 

51 k + 1 5 X K + 1 


( 68 ) 


Now the flows in (4 0 ) an( i (2) have identical directions, but the 
pressures and the magnitudes of the velocity are different. 

To assure equilibrium within (4), the pressures and velocity direc- 
tions in (4 0 ) and (4u) must be equal. And to satisfy this condition the 
assumed shock must be intensified by 

Obviously it shall 


= A0 f 


The new pressure in (4 0 ) 


is (according to section 2 b) 


(69) 


2 k 


Pi, = Ph ' + sin 2 7i A 7 j_ 


(70) 


where A 7 ^ is the change in shock angle 7 ^ and is computed by 
equation ( 55 ) 



(71) 


(For the calculation of C see section 2b) 

Posting pj^ = p^, equations ( 66 ), ( 68 ), (70), and (71) give 


kMc 


P 2 


1 - 


\j^2 c 


2B'\ = p. 


2k w 2 . 2 k - 1 


5U + 1 


Mj sin*^ 


K + 1 


p 2 k . „ A3 1 

p 3"3 n-r sin2r — 


(72) 


This equation is linear in A0 f . 
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Now the quantities 6' and. 7' can be defined 

6' = &i + £& ± = 6i + A9* (73) 


AO * 

7* = 7 ± + A7 i = 7i + (7*0 

C 


The pressure in ( 4 0 ) can be obtained directly from equation (70) . 

The Mach number Mij. itself can be determined according to chapter I, 

if 6' and y' are known. That in (4 U ) is likewise directly obtainable 
from M2 by the isentropic expansion A©' . 

The slight discrepancy between the values M^ o and is due, 

as stated in section J>, to the fact that the condition for pressure 
equality, owing to the change in static pressure after both shocks, does 
not require equal magnitude of velocity. So a small vortex layer along 
streamline FqS is to be expected. 

Before intersecting the expansion wave forms with the flow direction 
in zone (l) the angle (section 2) 


= + ^3 

2 


A0 


After intersecting the angle with the stream direction in zone (2) 
is 


Y = 


p 2 + - A©' 

2 


But there is a difference 5 between the flow directions in (l) 
and (2), so that the looked- for directional change is 


Pi + Px - A0 po + Pii - A©' 

£b = X + 5- Y= + & - 


2 


2 


(75) 
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The directional change of the shock front is 

^ c = 7 - (7' + 2©) (T^) 

A negative angle b and a positive angle c indicates that the 
expansion wave or shock front after crossing is in more downstream 
direction. 

For illustrative and comparative purposes, the graphical solution 
in figure 19 was made with the aid of the characteristics and shock polar . 
Here also the condition for pressure equality was replaced by velocity 
equality. 

The described mode of calculation is used in the following numerical 
example for illustration. 


Numerical Example 
The flow in zone (l) is: 

Pp/PO = 0.22905 Mp = 1.435 Up = 44.18P (See cascade 

example of the following chapter III, zone (3))- 

Data before crossing: 

intensity of expansion wave £& = 1° 

intensity of compression shock 6 = 3*03° 

shock angle of compression shock 7 = 47*91° 

With it the states in zone (3) become: 


p 3/P0p = 0.28478 = 1.469 


= 42.89° 


Therefore 

P2/P1 = 1-157 

p 2/ p 0p = 0.34560 M2 = 1-352 P2 = ^7 -75° 

Assumed shock intensity 

6p = 6 - A0= 3-03 - 1 = 2.03° 
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7± = ^5-29° 

P4o/p 0l - °-»55 


A - K * 1 M3 2 = 2.592 
2 


corresponding angle of shock 
pressure after shock 

Determination of the constants 


B = (lV^sinSyi - 1) = 0.092 


C = - sin^i 
Equality of pressure in (4 0 ) and. (Ay) gives 


sec^7| / — - 1] - sin 27 -M-|^ — 

1 U 1 1 B 2 


= 0.7504 


P 4 U = P 2 I 1 - 


KMo 

- A© l 


m 2 2 - 1 


t .,2 2 k . 0 A0 1 

“ P 4 0 - p 4 0 + p 5 m 5 rn Sin 27 i — 


which inserted gives 


A©' = = 0.89° 


Ay± = — i = 1.19° 
c 


After the crossing: 

shock intensity 6* = b^ + = 2.03 + O.89 = 2.92° 

y } = 7^ + A 7^ = 45«29 + 1.19 = 46.48P 


shock angle 
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Hence 

P4/P0 = 0.3303 n k = 47.08° 

Mi, = 1.3643 

M4 = 1.3640 

The comparison of the two Mach numbers indicates that the difference is 
quite small and lies within the calculation accuracy. 

Directional changes: 

Expansion Wave b = 43*036 + 3-034 - 47*47 = - 1*4° 

Shock front <) c = 47*91 - 46.48 + 1 = + 0.43° 


5* Crossing of Expansion Waves 

Each expansion wave is again replaced by n small waves. In fig- 
ure 20(a), two waves of intensity Z£)q and Z£) 2 cross each other in F. 

After crossing, the intensities are Z^B^' 311(1 20 2 ! . l n this case, only 

one stream direction is obtained in zone (4), when 


2©q + 202 = + Z®2 * (77) 

Application of the relations of section 2 results in 
p^ = p^ - Z^) 3 = p 2 - 2p 2 

that is 
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Since all other quantities are known, ZB 2 ' and A9i' can be com- 
puted from these two equations. 

The directional change of the Mach lines is like that in the pre- 
ceding section 


H-, + H* - AOq + Hk " A©/ 

b = -i | — + A 0 ^ - ± 1 ~ 

Hn + H-2 - ^©2 H5 + H4 " ^©2' 

c = + ZBo - 

2 2 

Since al 1 changes follow the same adiabatic curve, the condition for 
pressure equality yields equal velocity values at both sides of the 
streamline FS. Hence, no vortex layer will appear. Figure 20(b) 
represents the graphical solution. 


(79) 

(80) 


Airstream: 


Numerical Example 


p-^PO = 0.3295 M x = 1.366 

The intensity of the first expansion wave is: 


Hi =? 47.05° 

A©]_ = 0 . 99° • Therefore 


P2/P0 = °*3134 


M 2 = 1.402 


The second expansion wave intensity is: 
The conditions in zone (3) are then 


H 2 = ^5.50° 

A@p = 1.06° 


y Po = 0.31245 M 3 = 1.4041 


h 5 = 45.433° 

The two equations defining ZB ^ 1 and ASH,' are: 

AGD. + Z^ 2 = °-°3579 = + &2 


Kpjy 1 * 

P 4 = P 3 , ZBq = p 2 


\f*l 


KP2M2 

\J Mp 2 


zb 2 ' 
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hence 

£%' = 0.0184 = 1.054° 

AG£' = 0.0174 = O.997 0 
after which the conditions in zone (4) become: 

p^p 0 = 0.29722 M4 = 1.44 ^ = 44.01° 

By equation (79) and (80) the directional changes of the waves are 
b = + 0.5° and c = - 0.02° 


6. Reflection of Compression Shocks and Expansion Waves 

No difficulties occur in the determination of the conditions existing 
behind the reflected compression shock FB (fig. 21). Those in zone (2) 
can be defined according to chapter I, if the state of the airstream and 
the intensity 6 of shock AF are known. Obviously the reflected shock 
is of the same intensity as the impinging shock, so that the shock angle 7 
of the reflected shock and the conditions in zone (3) can be defined. 

The same holds true for the reflection of expansion waves, when the 
intensity of the expansion waves and their slope with respect to the wall 
are known. 
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CHAPTER III. THE CASCADE PROBLEM 


1. Problem 


Visualize a cascade of infinitely many and infinitely thin flat 
plates, of which two adjacent plates AB and A'B* are represented in 
figure 22. The angle of stagger is 90-0, the spacing t and the 
blade chord L. This cascade is exposed at angle of attack f to a 
supersonic flow Mp, pp. 

It is assumed that the flow is the same in all planes perpendicular 
to the plates and determines the force, that is, lift and drag as well 
as the pressure variation along the plate (blade) . 


To each plate there correspond interference lines (chapter I), that 
is, the expansion wave issuing from the leading edge and compression 
shock (fig. 22). 

At wide spacing, the separate blades of the cascade will not affect 
each other and the problem reduces to the single plate. 

Now if the spacing decreases for constant chord, the interference 
lines of one plate intersect those of the other, without, however, any 
force being exerted on the plates themselves for the time being. In 
this event, the force on each plate is the same as on the single plate, 
except that the wake flow is slightly disturbed. 

The values of t/L, below which the interference line of a plate 
begins to exert an effect on the adjacent one, are called (t/L) cr ^ 

"critical chord-spacing ratio." 


plates. After the crossings and reflections, new zones appear on both 
sides of the plate where the pressure as well as the velocities are 
unlike the uniform pressures and velocities to be found at either side 
of the plate. As a result, there is a change in the total force as well 
as the lift and drag on each plate. 

The mode of calculation consists in defining each intersection and 
reflection with the laws of chapter II and from it determining the con- 
ditions in the several zones. Integration of the various pressures on 
both sides of the plate gives then the total force, that is, the lift 
and drag. 


2. Method of Calculation 


At 



the interference lines are reflected on the 
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The resultant force still is perpendicular to the plate, hut no 
longer through the plate center, hence produces a moment with respect 
to the center. The position of the force is defined by statistical 
methods . 

This method is illustrated in the following example. 


3 . Example 

The cascade ABA'B' (fig- 23) with 30° angle of stagger, that is, 
p = 6o°, and at angle of attack of ijr = 3 ° is placed in a stream with 
M^ = 1.4004 (corresponding to v = 9 °) • 

The- blade spacing was assumed at the beginning, while the plate 
chord was so chosen after completion of the calculation that the expansion 
wave was reflected exactly once on the bottom side of the upper plate. 

It was found that t/L = 0 . 547 • 

The flow experiences a compression shock starting at the leading 
edge A'. The shock angle y - 49.57° is read from the shock tables 
and the shock front A 1 a can be plotted. 

Proceeding from the leading edge A, an expansion wave spreads out 
between the Mach lines Ax and Ay. The first forms with the airstream 
direction the angle = 45 . 56 °. The characteristics tables give 

My = 1.503, that is, the Mach number which is obtained at an expansion 

by 3° from the Mach number 1.4004. The corresponding Mach angle, that 
is, the angle which direction Ay forms with the plate, would be 
Py = 41.70°. Instead of the continuous expansion, assume an expansion 

in three stages, each corresponding to a 1° deflection. The conditions 
in zones 1 , 2 , 3 , and 4 are obtained from the characteristics tables, 
after which the directions Aa, Ab, and Ac can be defined. 

By applying the methods of chapter II to the calculation of the 
crossings a, b, c, e, f, g, l, m, n, p, q, s and the reflec- 
tions d, h, i, k, o, r, u, the static pressures, the Mach numbers 

(table 6 ), and the intensities of the expansion waves and compression 
shocks, as well as their directional changes (see table 7 and fig. 24) 
in the several zones, can be determined. 

The static pressures were referred to the standard stagnation pres- 
sure p 0l - 

The stagnation pressure changes were disregarded in the determina- 
tion of the Mach number. This change is rather small according to table 6 , 
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so that no appreciable advantage was to be gained by including it. (For 


The pressure distribution past the plate is obtained immediately and 
represented in figure 25 * There the passage of compression shocks and 
expansion waves is accompanied by a sudden pressure variation. Since the 
actual expansion is continuous, the serrated line is replaced by a smooth 
curve, such that the areas decisive for the force calculation are identical. 

Note that the pressures on both sides of the plate cancel out over a 
large portion of the chord. The resultant force can be determined by 
integration of the various pressures; the various spacings l are read 
directly from figure 23. 

The plate width was assumed at b = 1 . The result is 


calculation of 




= 0.2896 lower side 


= O.2963 upper side 


Downward resultant force 


K 


O.OO67 



lift coefficient 



drag coefficient 
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4. Calculation of Thrust, Tangential Force and Efficiency 

(a) The resultant force on the blade is resolved into two components. 
One - the thrust S - is normal to the plane of the cascade, the other - 
the tangential force T - parallel to it (fig. 26). 

If 


then 


K = resultant force per unit of area 
(90° - p) = cascade stagger angle 


S = K cos p 
T = K sin p 


(81) 


As functions of lift and drag 

S = A cos(p - i|f) - W sin(p - \|f) 
T = A sin(p - \|/) + W cos(p - \|f) 


Referring the force to the dynamic pressure of the airstream, 

gives the coefficients 




C G = COS P 

S 


c_ = — sin 3 
T qj. 


( 82 ) 


similar to the lift and drag coefficients, which can be obtained directly 
from c a and c w . 
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At fixed blade chord and fixed angle of attack the resultant force 
reaches its maximum value when adjacent blades do not affect each other, 
that is, when t/L > (t/L) crii -. this event 


where P2 1 is pressure at lower side (behind compression shock) and 
P2 is pressure at upper side (behind expansion wave) . 

For a given Mach number of flow and angle of attack the thrust and 
tangential force is maximum at p = 0° and p = 90°, respectively. 

At a given angle p and a given Mach number, T and S increase 
with increasing \|r. Owing to our assumptions may not exceed \|r s , in 

order to prevent subsonic flows on the bottom side of the plate. 

(b) Definition of efficiency (no friction): 

It is supposed that the air enters normal to the plane of the cascade 
at a speed v (fig- 26). The cascade moves with the tangential velocity 1 
and finds itself accordingly in a relative flow with an angle of attack \|r, 
whereby tan(p - \|/) = v/u. As a result of this flow, the two forces S 
and T normal and parallel to the plane of the cascade act on the plate; 

S and T are defined according to previous considerations. An efficiency 
is defined as on a propeller, by visualizing the blade being driven at 
speed u with respect to force T and so producing a force S in axial 
direction on the flowing air. Then the power input is T X u, the power 
output S x v and the efficiency is 


The efficiency is seen to be dependent on \|f and p only. At 
constant p it decreases with increasing \|r. At \|r = constant, t] has 
a maximum, if 


K = p 2 ' - p 2 



(83) 


or 



tan p 1 + tan ^ tan p 


^L= 0 
8(3 


(84) 
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that is, when 


tan p = tan + \/ ( tan + 1 


\/( J 


which approximately gives 


p = 45 0 + - ♦ 
2 


( 85 ) 


The maximum efficiency is then 


^max 

\|/ = Constant 


1 - tan — 


1 + tan — 
2 


( 86 ) 


At small values of >|r, tan \[r = \|r and \|r^ is negligibly small, hence 


at p = 4-5° + — + 
2 


‘max 


1 - » 
1 + ♦ 


( 87 ) 


The efficiency for various p and \|r is represented in table 8 and 
figure 27 '- 
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CHAPTER IV. LINEARIZED CASCADE THEORY 
1. Assumptions 


The theory is based upon the following: 

(a) All disturbances are small in the sense that all interference 
lines may be regarded as Mach lines. The expansions are simply concen- 
trated in a Mach line and the compression shocks replaced by Mach com- 
pression waves. 

(b) Intensity and direction of waves are not changed by intersection 
of expansion and compression waves. The justification of this assumption 
is indicated in the preceding numerical example, where it was shown that 
the directional changes of the wave fronts are small, as a rule. 

On these premises, the interference lines AA 1 and AA" parallel 
to BE 1 and BB" start from the small disturbances A and B 
(fig. 28(a)). At the intersection in a the directions of the waves AA 1 
and BB' as well as their intensities remain unchanged. The pressure 
and the velocities in the zones (2), (3)> a 11 "! (4) are defined by the laws 
of chapter II. In the holograph these assumptions imply that the char- 
acteristics network in the applied zone is replaced by a parallelogram 
(fig. 28(b)). 


2. Linearization of Cascade Problem 

The application of these simplifications to the solution of the 
cascade problem produces parallel Mach lines within the cascade, which 
remain parallel after crossings or reflections (fig. 29(a)). (L = plate 
chord, t = spacing and \|» = angle of attack.) The Mach lines Aa 
and A' a emanate from the leading edges A and A'; the angles aA'X' 
and aAX are Mach angles and both equal to On passing through A' a, 

the flow experiences a compression and a directional change \|r, along Aa 
an expansion with the same directional change. 

The pressure in (2) and ( 3 ) can be defined by the laws of isentropic 
expansion and compression (chapter I); that of zone (4) is computed the 
same way from the pressure in (2) and is obviously equal to p 1 , as seen 

in the holograph (fig. 29(b)). But the flow direction in (4) differs 
from that in (l) by an angle 2i|r. 

The Mach line aC' intersects the plate at C* and is reflected 
along C'E, whereby C'E is parallel to A' a. The pressure in zone ( 5 ) 
is again equal to that in (3) and the flow is obviously parallel again 
to the plate. 
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On passing through DE' the flow from (4) and (6) is compressed - 
the reflected wave DE' - so that in (7) the direction and the velocity 
of flow are the same as in (l)$ the same applies to the flows in (6) 
and (2) . 

Thus it is seen that the corresponding zones repeat themselves, hence 
that the further conditions are completely known without new calculations. 


3 . Calculation of Lift and Drag 

The pressure variation on either side of the plate can be plotted 
(figs. 29(c) and 29(d)). The pressure remains constant over the lengths 
AC, CD, DE, EF and FB and over A'C', C'D', D'E', E'F' and F'B' 

where the interference lines strike the plate. 

Along CD the pressures on both sides are equal and cancel out, 
whereas a downward pressure difference p^ - pp, obviously perpendicular 

to the plate, acts on AC and EF, and an identical upward pressure dif- 
ference on DE . 

The pressure pattern in figure 29(e) repeats itself in length direc- 
tion of the plate over the period L^. If the plate chord is chosen 

exactly like or a multiple of it, there is no resultant force, that 

is, a plate of this length has neither lift nor wave resistance. 

For the values of L, which satisfy the inequality 
L 0 < (L - nip) < (Lp - Lo) 

whereby n can be = 0, 1, 2, . . ., the resultant force reaches its 
maximum value, and then 

K = (p 5 - p 2 )Io ( 88 ) 


Hence it serves no useful purpose to make the plate longer than Lq, 
because there is no more lift increase anyhow. On the other hand, a 
moment occurs and, in the presence of friction, the drag would increase 
unnecessarily. The boundary Lo(= AC) is the plate length not touched 
by interference lines of the other plate and can be defined geometrically 
in terms of cascade spacing t and angles 0, \jf, and p. 
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L 0 = t 


sin[p - (up + Y)J 
sin(^-i_ + \lr) 


= 2Lq + t 


sinjp + (^ + \|f)] 
sinCn! - \|f) 


(89) 


Accordingly the best ratio of spacing/chord is 


t t 


sin(p 1 + i|r) 


L LO sin^p - (|i-^ + \|<J 


(90) 


Now c a and can be determined when 


1. L = nL^ then c a = c w = 0 


2 . Lq < (L - nLi) < (Lp - Lq), "the boundary values are 


K 

C a = — cos I|r = 

dl L 


P 5 - P 2 


F cos \|r 


K P3 - P2 

Cy = sin = — F sxn i|r 


(91) 


where 


F = — 


sin[p - ( Pp + ^)] 
L sinCpj^ + \l») 


3- L 0 > (L - nLp) 


c a = 


C V “ 


p 3 ~ p 2 
P 3 - P 2 


{L - nLi)cos \|r 


q.]L 


(L - nLi)sin \|r 


( 92 ) 
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4. (L - nli) > (Ll - Lo) 

•J^(L - nil) - (Li - Lq)Jcos \|r 


Co — 


p 3 - p 2| 

q x L 


~w 


Px - Pp r 1 

= -2. £ (L - nil) - (Ll - Lo) sin 

q x L L J 


(93) 


The linearization can be extended to the pressures P2 and p^; 
admittedly then only when the angle of attack is sufficiently small. 5 

The pressures can be defined by the laws of small variations 
(chapter II). Thus 


P2(3) = Pi 1 f 


kM-, 2 ^ 


V 


(94) 


M ± 2 - 1 


Inserting these values in the above formulas for c a and c w , while 
expressing the dynamic pressure with 

91 = | k Pi Mi 2 

and the values 1 and i|f for cos i|r and sin ty, gives as for the 
isolated plate, 

3ln the following table the pressures after expansion of p 0 = 0.31404 
(corresponding to Mp = l-4004) are represented in terms of the expansion 
angle: 

pp = pressure according to isentropic law of expansion 
p = pressure according to the laws of small variations (chapter II) 

t—±j 

Mr P n ^ _ 


p 2/P0 

P2l/p0 

( p 2 - p ?T. ) /Pp p ercent 0.1 


0.29906 0.28478 0.27114 0.25809 0.23363 
O.29865 0.28335 O.26718 0.25266 0.22196 


0.5 


1-5 


2.1 


5.0 
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The factor F approaches 1 when t/L = t/Lg. The theory is now illus- 
trated on the following numerical example. 

4 . Numerical Example 

The cascade of the numerical example in chapter III is applied again 
with the same airstream as by linearized theory, figure 30. 

It was 


t/L = 0.547 

p = 6o° 

M x = 1.4004 

II 

/ Po = 0 . 3i4o4 

m = 45 . 56 ° 


The Mach lines within the cascade can now be plotted. By equation ( 89 ) 


L 0 = O.lVkL 


Geometrically defined are 


(L x - L 0 ) = O. 778 L 


so that 


(L - Li) = 0.078L 
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The tables of characteristics give 


p 2 /p 0 i = 0.2711 (expansion by 3° starting from P 1 /p 0l ) 
P3/P0 = 0.3640 (isentropic compression by 3 °) 


Assuming the plate width at one cm, gives: 


resultant force 


K 


= (L - Li) — — = 0.3436 

p 0l p 0l 

resultant force per unit length = O.OO67 


lift coefficient c a = O.OI56 
drag coefficient = 0.0008 


The pressure distribution on both sides and the resultant pressure 
are shown in figure 31- 

5. Comparison With Exact Method 

Instead of the lengthy calculations of all crossings and reflections, 
the linearized theory affords a quick and simple solution of the cascade 
problem. At small angles the results are reliable and the errors small, 
as seen from the comparison with the numerical examples in section 3, 
chapter III and the preceding section. 


c a( exact) - c a ( linearized) _ _ 2 percent 
c a( exact) 

The interference lines of the linearized solution within the cascade 
the Mach lines - are included in figure 23 for comparison. It is seen 
that the zones governing the resultant pressure are smaller by linearized 
theory. 

The pressure distribution of the linearized example is also shown 
in figure 25- 
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CHAPTER V. SCHLIEREN PHOTOGRAPHS OF CASCADE FLOW 
1. Cascade Geometry 


A disturbance in supersonic flow is known to spread out only down- 
stream of the source of disturbance. So the pressures and velocities 
on one of the sides of a profile, stipulated by the form of the surface, 
are not influenced by the other side. 

This property is used to represent the flow through a cascade con- 
sisting of a number of infinitely thin plates. Two profiles with a flat 
surface on one side are so assembled that their flat sides face each 
other and are parallel. The flow between the parallel sides is exactly 
the same as that between two adjacent plates of the cascades. 

The two profiles can be moved apart or shifted relative to one 
another, so that any desired ratio t/L and any stagger angle can be 
obtained. 


The experimental cascade was patterned after the cascade in the 
numerical example of chapter III, which had the same angle of stagger 
of 50°. The Mach number of flow was - as in previous calculations - 
M = 1.40; the spacing ratio was t/L = 0.517- The angle of attack \|r 
ranged from 0°, l-5°> 3° 1° 4-5°. 


The maximum profile thickness was so chosen that no blocking of the 
tunnel (section 2) was produced at the selected Mach number and that the 
deflection of the profiles at maximum angle of attack is small. 

Now at M = 1.40 the deflection due to compression shock, which 
exactly leads to sonic velocity, is 5 S = 9°- As there is to be no sub- 
sonic flow in the test section and since the angle of attack was assumed 
at 4*5°, the leading edge of the profile may at most form an angle of 
about 4°, which corresponds to the constructed profile. 


The compression shock is not separated at the leading edge of an 
infinitely thin plate or an infinitely sharp wedge of sufficiently small 
included angle. Therefore the leading edge shall be as sharp as possible. 
It succeeded in attaining a thickness of 0.05 — 0.07 rani 4 " so that the 

distance of the separated shock from the edge is scarcely visible. 


The profile chord L was 118 mm, so that the cascade lies within 
the tunnel window. Since the tunnel itself was 400 mm wide, the width 
of the profile was limited to 398 mm, figure 32 . 


k ” 

ulurit Co., Affoltern, Zurich. 
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2 . Experimental Setup 

The previously described profiles were mounted in the test section 
of the supersonic tunnel of the Institute5 on four supports (fig. 33)* 

The compression shock issuing from the leading edge of the top 
profile could not be reflected at the upper tunnel wall at maximum \|/, 
because the deflection to be made retrogressive at the wall was too great 
for the Mach number prevailing behind the shock. To avoid blocking in 
this region, a bend had to be made in the upper nozzle wall (fig. 3*0 • 

The position of the bend was so chosen that the fan of expansions emanating 
from it hits the cascade downstream from the entering edge. This adjusts 
the wall to the flow direction after the shock to some extent as well as 
raises the Mach number between the upper plate and the nozzle wall. 

The Mach number in the test section before the cascade was deter- 
mined by pressure measurements at the upper, lateral, and lower walls. 

The investigation was carried out at a moisture content of air of about 
0.5 g water/kg air. 


r 

3. Schlieren Photographs 0 

The schlieren photographs illustrating the flow through the plate 
cascade at \|r = 0°, 1,5° and 3° are represented in figures 36, 37, 
and 38. Since a conical jet regime is involved, the photographs appear 
as shadows of the profiles. Figure 35 shows the position of the optical 
axis with respect to the cascade; it is seen that the shadows of the pro- 
files are distorted on the mirror. At the top profile the perspective 
effect is more obvious, because the optical axis is closer to the bottom 
profile . 

The equality of Mach f s angle in figure 37 U = 0 °) is indicative 
of an unchanged Mach number in the cascade. The visible disturbances 
within the cascade may be due to the fact that the plate surfaces do not 
exactly agree with the flow direction, or to thickening of the leading 
edges by a boundary layer. 

In figure 38 (\|r = 1 . 5 °) the interference lines inside the cascade 
are almost parallel, as stipulated by the linearized theory. 


5 See Report No. 8 of the Institute for Aerodynamics, at the E.T.H, 
ref. 1. 

^For description of schlieren apparatus see Report No. 8 of E.T.H. 
Institute . 
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At i|r = 3° (fig* 39) the deflection of the shock front at crossing 
of .the expansion wave emanating from the top leading edge is plainly 
visible. Figure 4-0 represents an enlargement of the crossing to illus- 
trate the numerical example in chapter III. The interference lines inside 
the cascade for this example are again shown in figure 4-1 at smaller scale 
(compare also fig. 23 ), whereby the perspective effect is indicated. 

In the majority of photographs the retardation of the flow near the 
tunnel wall leads to separation of the head waves. 

The flow in all photographs is from left to right. 
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CHAPTER VI. THE FLAT PLATE CASCADE AT SUDDEN 
ANGLE-OF-ATTACK CHANGE 
1. Problem 


Visualize a cascade of flat plates in a flow with relative velocity W 
at an angle of attack ty. A supersonic flow which may be regarded as two- 
dimensional prevails throughout the cascade. At a given moment the angle 
of attack of the airstream changes from to V within an infinitely 
short time interval. The transition to the new state , which is to last 
for a period, is analyzed. 

Such a change in the angle of attack takes place when the cascade 
moves in an absolute flow which has not the same speed at every point, 
or when one of the velocity components of the flow, normal or parallel 
to the plane of the cascade, varies with respect to time. 

Resolving the velocity W in two components V and U (fig* ^-2) 
normal and parallel to the plates, the change of the angle of attack, 
small in itself, can be regarded as a change of component V. This change 
in V is obtained by superposition of a velocity vq, which has the same 

direction as V and is obviously small compared to V and consequently 
smaller than sonic velocity. From the assumption of a small angle of 
attack, it follows that velocity component U remains greater than sonic 
velocity. Besides, an eventual variation of this component U is 
disregarded. 

The problem therefore reduces to the study of the new forces on the 
cascade, resulting from a gust v 0 7 which, together with the velocity U 

enters perpendicular to the plates. 

Biot (ref. 5) solved the problem of an isolated plate by means of 
"unsteady sources." This method is applied to the cascade problem. But 
first the unsteady source is described in more detail. Since the plates 
are to be partly replaced by such sources, the pressures and velocities 
originating from a source distribution are analyzed. Then Biot f s results 
for the isolated plate are correlated and extended to the cascade. The 
special case of straight cascade (nonstaggered) is examined. 


7 By "gust" is meant a continued, uniform vertical velocity distri- 
bution vq. 
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2. The Unsteady Source 

According to linearized theory, the general potential equation (5) 
for two-dimensional unsteady flow can be simplified to 


gap . m 2 ) + - a SI .1^,0 (96) 

3x 2 3y 2 a Sx 3t a2 3t 2 


cp = flow potential. 

For a system of coordinates moving with velocity u(u/a = M), that 
is, air at rest at infinity, this equation gives the acoustic wave equa- 
tion for two-dimensional motion 


+ at _ ^ at = 0 (97) 

3x 2 3y 2 a 2 3t 2 


One solution for a linear sound source is 


cp = k cosh" 1 ^ (98) 


where r = 
velocity. 

This 


\Jx^ + y^ and K = constant with dimensional length times 


solution is rewritten in the form 



(99) 


= -k log. 




1 at + \/a 2 t 2 - r 2 


J 
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It represents a cylindrical wave varying in time rate. At t = r/a, 
cp = 0, that is, if such a singularity appears in the zero point of the 
coordinate system, its effect is diffused inside a circle of radius 
r = a-t. 

If such a source appears at the point (x,0) - on the x-axis - at 
period tp, the potential in a point P(x,y) of the surroundings of this 

source at a given period (fig. 43) is 


(100) 

In this case 

r = ^(x - xi) 2 + y 2 

and the following velocity components are obtained by simple differentiation 


cp = k log e 


a(t - tp) + \J a^(t - tp) 


2 2 


Sep ,1 a(t - t^) 
v r = — = k 

dr r yk 2(t - ti)2 - r 2 

_ &£_ = k ( x - x l) a(t - tx) 

^ x r2 ^ a 2(t - ti)2 - r 2 


(101a) 


(101b) 


v y = ^ = k x a(t ~ tl) — ( 101 c) 

^ r \j^( t - ti)2 - r2 

When y is small compared to a(t - tj_) - near the source - the for- 
mula (lOla) becomes 


v r = k/r 
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the same as that of a steady wave in incompressible flow, hence with Q 
denoting the strength of the source (dimensional length x speed) 


= ± 
2 * 


The pressure in the same point is computed by 


= 


<dp _ paQ 1 

2:rt \ja 2 (t - tp) 2 - r 2 


( 102 ) 


It will be noted that vy always equals zero for y = 0, except at r = 0 

in the source itself. It means that such a source delivers at no other 
place on the x-axis a velocity component parallel to the y-axis. 


3 . Pressure and Velocity of a Periodically 
Arising Source Distribution 

Consider a continuous distribution of infinitely small sources over 
the length 0A (fig. 44) along the negative x-axis. The distance 0A 
increases linearly with the time: 0A = Ut, where U is a constant 

velocity and the sources on the x-axis appear momentarily at the point 
where A arrives at the moment. The strength of this source distri- 
bution per unit length of 0A is assumed equal to q (dimension of a 
velocity) and remains constant in time. 


(a) Pressure 

At point P(x,y) (fig- 44) the pressure p of the source distri- 
bution at time T is, by equation (102) 


P = 


paq 


>xp=0 


dXn 


27t X 1 - tp) 2 - (x - xp) 2 - y 2 


( 103 ) 


ty the time of origin of the source in point xy. 
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With the following variable transformation 


i = 


X 

at 





£l = 


X 1 

at]_ 


a . 1 

- = sin u = — 

U M 


(104) 


we get 


P 

paq. 


1_ 

2jt 


-l/ £ 




Sin ^ \J( 1 + Cp sin p ) 2 - (i - Cp ) 2 - rf 


(105) 


The boundaries should be defined before the integral is evaluated. 
For the function in the denominator is real only in the zone affected 
by the source distribution; this is bound by the Mach line AM and the 
circle with center 0 and radius a-t. Hence the integration must be 
made between the zero places of the function where it is real. 

Posting 


(l + Cl sin p ) 2 - (C - Cl)" - - n 2 = 0 ( 106 ) 


the new boundaries are found at 


Si 

Si 


(1) 

( 2 ) 


> = 


-(c + Sin p) - \j( 1 + sin p) 2 - tj 2 cos^ 


(107) 


-cos2p 


To 

tion of 
results 

Q(S,Si) 


get an idea of the integrating process as function of the posi- 

point P, Ci( 1 ) and C^ 2 ) are plotted in terms of £. It 

in two curves of the second degree, which cross in point 
(fig. 45(a)) whereby 
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i = - 



T] COS |I 

sin \i s 


Si- 


' T) - COS 
ysin |jl cos 



(108) 


The shaded area represents the r^nge in which the integration should 
he made. At small £ values up to £ = ^1 - r^ 2 , integrate between 

and and- then between and the £ axis. In figure 45 (b) 

the integral limits are shown p ,ted in the x,y-plane for explanation. 
The reason for not integrating over positive values is the absence 

of sources in the right-hand half plane. 

Two integration cases are differentiated 



and 


- \/ a^t^ - < x < + \ja£t^ - y^ 

y 


In the first instance the pressure integral is 



(no) 


paq 2jt 
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But as 

the root, it 


and £-^(2) are the solutions of the expression below 
can be rewritten as 


P 1 r ? i (2) d? i 

” ' ■ /(,»' - - fc) 

With the substitution 

6l - £i (l) + ^i (2) - 5i (l ^sin 2 9 


( 111 ) 


( 112 ) 


this integral 


8 


gives 



paq 2 cos p 


a formula that is independent of t) = y/at. 

In the second case, if point P is so situated that 



it results in 


P 

paq 





/(l + ^ sin n) 2 - (^ - ^) 2 - t \ 2 


(115) 


(11A) 


®As long as the function in the denominator can be brought, with 
the aid of the integral limits, into the form of equation (ill), the 
integral gives the same value. 
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which after evaluation gives 


p y _ £ 

paq at cos p 


cos 


-1 


(S + sin p) 



+ £ sin p) 


2 


2 2 
T) COS p 


(115) 


The ensuing pressure pattern along a line y = Constant is repre- 
sented in figure 46. For each y the pattern consists of two pieces. 

In the first piece the pressure is constant and equal to p c , that is, 

along the length EF between the points where the Mach line emanating 
from A and the circle with center 0 and radius a-t intersects the 
line y = Constant. The second piece is composed of length 

FO' = - \j a^t2 - y2 and O’G = + \/a2t2 - y2, where the pressure is vari- 
able; at G the pressure is zero. At y = y c the constant portion 

disappears and wherever y = at, the pressure becomes zero. 

At y = 0 it represents Biot’s case with the integral limits 


that is 


1 

sin [i 


< i 


< 1 


at 

sin \i 


< x 


< - at 


and -1 < £ < + 1 


> 

and -at < x < + at 


(116) 


The pressure p c has the same value as before 


£c_ _ 1 

P a Q. 2 cos p 


but the second piece of the pressure distribution becomes 


p y=o 

paq 


/ 


2rt cos p 


cos 


-1 


at 


+ sin p 


1 + — sin p 
at 




(117 ) 9 


9q corresponds to Biot's 2 Vq. 


(118) 


NACA TM 1369 


57 


It should be noted that a pressure effect appears also outside the 
area in which the sources are distributed, because the source in 0 
affects the area inside the circle a-t as mentioned before. 

(b) Velocity Calculation 

In general, the velocity component is defined by the integration 
of the portions stemming from a single source (eq. (101)). In our problem 
the velocity Vy is of particular interest. 

It becomes 


1 

Vy 2* 


r x r=o 

'xi=-Ut 


a(t - tp)y dxp 


(x - X-.T + t 


a 2 (t - t 1 ) 2 




x ,) 2 + f 


If r = \J(x - xp) 2 + y 2 is small compared to a(t - tp) , that is, 
for the places close to the x-axis, this equation simplifies to 



Letting y approach 0, positive y, the results for negative values 
of x are 



( 120 ) 


which may "be designated by Vq (as in Biot's report). 

For positive x values, Vy = 0. 

It indicates that such a source distribution gives a uniform vertical 
velocity Vq over the distance of the x-axis where the sources are. (For 
negative y, inverse velocities result.) 
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Biot mentioned this fact in his report and used it to calculate the 
pressure distribution over a plate in a vertical gust (compare next sec- 
tion) . 

The same variable change as in the pressure calculation gives 

( 2 ) 


v 0 

v v = — 
y n 


rtf 

SiU) 


r] ( 1 + Cq sin p)d£q 


[(C - Cq ) 2 + n 2 J j^(l + Cq sin p) 2 - - C x ) 2 + T] 2 ] 


( 121 ) 


The arguments for the integral limit are the same as for the pres- 
sure integral and Cq^"^; Ci^ 2 ^ is given by equation (107). Integrating 


between Cq an< i Cq^ 2 \ that is, when 


'l - r) cos p' 
sin p 


< i < 




it is seen that the integral gives the value — , so that v v has the 

r\ y 

constant Vq for this range of £• 

For the second case \[i - n 2 < C < +\J 1 - T)^ ^ the integration 
between Cq^ and 0 indicates that^ 



( 

_ 



at - tan~^ 

ii\yi - (c 2 + 7 2 ) 


at 

sin p(C 2 + T) 2 ) - C_ 







( 122 ) 


This identifies the velocity distribution on the lines y = Constant 

(fig- Vf). 


10 , 


See note on p. 66. 
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From the calculation of the pressure and velocity distribution over 
the lines y = Constant, it is apparent that the flow outside the circle 
of radius a-t is steady. This is true from the physical standpoint too, 
since the gust front does not affect this area. 

4. Single Flat Plate in a Vertical Gust 
(Biot 194-5) 

The flat plate AB of length l at supersonic velocity U enters 
a gust with the transverse velocity v Q (fig- 48(a)). Since the trans- 
verse velocity on the plate must be zero (no flow through plate) an equal 
and opposite velocity is superposed on the gust velocity vq in place 

of the plate. This velocity can be visualized as reflection of the gust 
on the plate (fig. 48(b)). 

Since velocity U is greater than the sonic velocity, the sides of 
the plates are not affected by one another, so that one side of the plate 
can be analyzed separately. The pressure acting on one side is exactly 
the same as on the other, except with inverse prefix. As the interference 
velocity Vq is much smaller than velocity U, the linearized potential 

equation can be applied to the stream potential. Selecting a system of 
coordinates that moves with the velocity U, (eq. (97)) according to 
which the disturbances are diffused with sonic velocity, can be applied 
to the flow potential. 

Biot's method replaces the part AO of the plate struck by the 
gust, by unsteady sources. This source distribution, which increases 
in time, yields a uniform velocity Vq normal to the plate, hence sat- 
isfies the boundary condition on the bottom side of the plate. 

If the plate enters the gust at time t = 0, the distance at time t 
is AO = -Ut, the origin of the coordinate system being located in the 
gust front. 

The results of section 3 can be applied directly, and the pressure 
variation along the plate defined (fig- 49(a)). As it is dependent 
solely on x/at the patterns are like those for the different Mach num- 
bers. The total force on the plate - the lift - is obtained by integra- 
tion of the pressure pattern. Three phases are involved here (fig- 49 ) • 

I (U + a)t = l 

that is, the trailing edge is outside the effective range of the gust 
front; 
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II (U + a)t> l > (U - a)t 

that is, the trailing edge is inside the effective range of the gust 
front j 

III (U - a)t = l 

that is, the entire plate is outside the effective range of the gust 
front, hence is no longer exposed to any unsteady effect. 

The integration gives the following lift values of the three phases: 




Al 

ut . 

a 

4 - 



(125) 



2pav 0 t 

l sin 

M- 



A II 

1 -1 

1 (l ■ 

ut 

c 



sin'l 

1 (~ 

\ 

_ i 1 + 5. 

2pavQ l 

Jt COS [1 

sin 

V-WiCJ M- 

l ) 

+ 

iti 

sin p\Ut 

J 2 


(124) 


^The integral 


U - a) \/ i - l 


1 = 


appears in the calculation of Aj and Ajj. With no boundary, the solu- 
tion is 


T • -±y A . ,/l + A£ 

I = sin -% - sin-l 1 

\/a 2 - 1 \ A + £ 


which gives 


I = it 


between the limits -1 and +1. 


\ 


\l A 2 - V 


ro 1 ^ 
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sm 


-1 


to be taken between 



and 


♦i- 


a iii 1 

2 pavQZ cos p. 


(125) 


In phase I and II the lift increases continuously with the time and 
reaches a maximum in phase III, where it becomes independent of the time. 
In the last phase the lift is the same as on a plate at angle vo/u in 
steady flow. 


5 . The Straight Cascade 

The cascade problem is unlike that of the plate to the extent that 
the plates mutually interfere. The sources replacing the portion of the 
plate struck by the gust create a pressure on the adjacent plates. They 
also produce a velocity v u , which in order to satisfy the boundary con- 
dition of no through flow of the plate, makes a change in that source 
distribution necessary. 


Since the disturbances are small the solution of the single plate 
can be superposed in the sense of the linearized theory of the adjacent 
plate effect. 

As shown in sections 2 and 3, the unsteady source - and the source 
distribution - which lies on the x-axis, produces no vertical velocity 
component along this axis, outside the distance, where it is. 

This characteristic enables the velocity component v y to be 

replaced by an additive source distribution along the particular parts 
of the plate, which gives the velocity at each point. The new sources 
create a further pressure on the plate itself - and in general react on 
the adjacent plates. 

The total force - the lift - on each plate consists then of the 
lift of the undisturbed plate (A u ), the lift from the pressure p y of 

the sources of the adjacent plates (A&) and the lift (A v ) of the new 
source distribution due to velocity v y . 

Suppose that h is the plate spacing and L the plate chord of 
the straight cascade ff (fig. 50). The lines AM . . . represent the 
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Mach lines emanating from the leading edge, where angle MAB is the 
Mach angle p = sin - "*" u/a. 

Then the following approximation is made: the relative flow and 

the plate form in reality the angle = tan“"K Vq/u) . But as vq is 

small compared to U, this is negligibly small with respect to p, 
and it can be assumed that the Mach line itself rather than the relative 
flow direction forms the angle p. In this event, the Mach lines form 
the same angle with both sides of the plate. The Mach lines emerging 
from the leading edge strike both sides of the plate at the same distance 

AE from the leading edge. At (h/Z) > tan p the points are not located 

on the plates, and the plates do not influence each other. Consequently 
the cases where (h/Z) < tan p are examined. 

At time t = 0 the cascade is directly in front of the gust; the 
origin of the coordinates is placed in the gust front. In the first 
time intervals of the phenomenon the disturbances have not spread out 
enough to be able to influence the adjacent plates. As in figure 51 , 
the distance is a-t < h, so that the circles with center 0 and 
radius a-t do not touch the plates. Lift and pressure distribution 
are the same as on the single plate. 

As soon as t > h/a, the plate AB comes within the effective range 

of its adjacent plates. On EFG (fig. 52) the source distributions A ! 0 ! 
and A"o" create an additional pressure which can be computed according 
to section 3* The points F and G are then the points of intersection 
of both circles with center 0 f and 0” and radius a-t with plate AB. 
It is readily apparent that the additive pressure on EF is constant and, 
according to equation (113), has the value 


Pc _ 1_ 

pavo cos p 


(126) 


If t] = h/at is inserted (y = h) in equation (115), the pressure 
on FG follows at 


Ph 1 1 

1 

(St + sin t 

1 

pavQ rt cos p 

/(‘•ft 

sin ^ 

2 ^2 2 

1 -a^ C ° B * 


The ensuing additive pressure is represented in figure 52- 


( 127 ) 
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In addition, the following condition must be satisfied: The normal 

velocities Vy = h created by the source distribution A' O' and A o' 

are reflected on EO, so that at that point the gust is partly compen- 
sated. The source distribution to be applied is to compensate the veloc- 
ity (vq - Vy) . The velocity Vy is computed as in section 3> and the 

pressure py - along the particular plate - is obtained by integration 
of the pressure contribution of each source. Assuming the local veloc- 
ity Vy on a small distance dXy to be constant, the yield of the source 

distribution per unit length on this small distance is then q = 2v y . 

Along this area of the plate the source distribution produces the pres- 
sure (compare section 2) 


4Py=o = 


paVy 


cLx-i 


\/a2(t - t x ) 2 - r2 


( 128 ) 


hence 


P 


v y=o 



pa v y dxi 

\^a2(t - ti) 2 - r2 


(129) 


with Vy periodically and locally variable. It is best to solve the 

integral graphically for each particular case. The arguments for the 
integral limit are the same as before. The lift contribution Ay at 

any instant is obtained by integration of the ensuing pressure plot. 

To obtain the resultant pressure, this pressure is superimposed on 
the two previous pressure distributions. 

In the following, the pressure contribution due to the additive 
pressure p^ is calculated. Three phases, depending on time and 
ratio h/Z, are involved: 
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I. When l > 


I^Ut + ^a^t 2 - h 2 ^ (fig. 52 ), the lift is 


p-\ya 2 t 2 -h 2 ~ + \/a 2 t 2 -h 2 

Ayi = 2 // v \ Pc dx + 2 / / P h dx 

yA J (-Ut + ^— j J a 2 t 2_ h 2 n 


with the previously employed variable change and with y = h 


V ~ 


r - n/i-b 2 

2at / 

J _ Ut+y/tan 

at 


i£ 


p c + 2at 



P d£ 


■N^ir 




= 2atp c (£) , + 2at(p£) 

C . Ut+y/tan p _ x/II^2 

at 


at 
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By equation (127) 


d P h 


d£ 


P av 0 (l + £ sin p) - T] 2 

£(l + £ sin p) 2 - T]2 C os2pJ \Jl - £ 2 - T| 2 


(131) 


hence the integral 


I = 


.. dp h 


-x/l^r 


£ — d£ 
d£ 


^/l^i 2 


(l + £ sin p - r|2)£ 


- \j 1-T) 2 \^i - £ 2 _ t]2 (1 + ^ sin p) 2 - T} 2 cos 2 pJ 


d£ 
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Its evaluation gives 


I = - 


it 


sm \i cos 1 1 


(COS [i - 1) 


x + Ti(cos 4+1) 

T) + COS |! 


(132) 


consequently 


Vi 


Spa^UQ 


2^1 sin p cos p 


f 


1 - T) COS p) + (COS p - 1) 


1 + 


T]( COS p + 1) 


(cos p + q) 


II. But if (fig. 53) 


(133) 


^Ut + \ v/ /a2t2 + y2^ > l > (ift - 


a 2 t 2 _ y2 


the evaluation of the integral gives the formula 


Ayii ~ 2 Pe I 


l - Ut \ 2pavp 

at j cos p 


'Z 


1 - if 


2pavp 


tc sin 11 cos p. 


sm 


•% - 


sin-1 + iA . _ _A_\ 

\/a 2 - 1 A+S 2^ ^A 2 - 1 1 


(tj - cos p) 


sin"% - 


B 


V^B 2 " 


• -i 1 + BS jt /-, 
sin x + —/l - 

B + S 2 ' 


B 


\J & 


(134) 
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where 


S = 




l - Ut 

a2t2 _ y2 


A = 


1 + T] COS p 
sin |jl \/l - t)2 


and B = 


1 - rj cos p 
sin | i\J 1 - r|2 


The pressure p e is obtained from equation (127)> when x 
is inserted, at 


Pe 


cos 


-1 


pavQ tc cos p 


l - Ut 
at 


\ 


+ sin p 


1 + >---!£ sin p\ 2 - -Si- oos2,i 


V 


at 




III. If l ~ (U - at), py = p c along the entire distance 
that the additive lift 


^ylll = 1 L _ h N 

2pav 0 cos pi tan p 


reaches a value that is independent of the time. 


Note on the Velocity Integral 

By a simple transformation the integral can be rewritten in 
following form: 


I = 


■Si 


( 2) ( A + B £l)<% 


J t- 


(q) (£q 2 + + B£q + y 


( l - Ut) 


(135) 


EB, so 


(136) 


the 


NACA TM 1369 


67 


This integral can he solved by means of tables (integral table, 
Part I, by Grobner and Hofreiter, Springer Co., Wien). Although the 
general solution is quite complicated, the result is found to be inde- 
pendent of once the limits have been inserted. 


Bearing in mind that the integration limits are the solutions of 
the function below the root, the integral is rewritten as function of 
the limits 


X-/, 


Cl 


( 2 ) 


(1 - £1 sin 


Cl 


( 1 ) 


(1 - 5 1 sin td 2 - cos 2 ^ - C 1 (l) )(c i (2) - tj l/(Ci - Ci (l) )(Ci (2) - Ci) cos 


The following substitutions are made consecutively: 


1. X = 1 - £y sin (i 

The limits are thus Xy and X2 (X 2 > Xy) 


2. Y 


y JXyXg_ 

X Xy + X 2 


1 


3. Z = Y(X 2 + Xy)/(X 2 - Xy) 

k. t = z 2 

4X 1 X« - cos2|i(x n + Xp) 2 
5. s = — — A' = — — — — — 

A' + t cos 2 p(Xy - x 2 ) 2 


6. Numerical Example 

Dimensions of cascade h/l = 0-55 

Mach number of flow M = l.klk (= v/2) corresponding to a 

H = ^5° 


Mach angle 
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The period t e 


up to the end of the phenomenon is determined by 


(U - a) t e = l 


whence 


"^e 


l_ 

0.4l4a 


All time intervals are referred to y/a in order to obtain a dimen- 
sionless ratio (= 1 /t)) . Then 


- 7 - = ^-38 (y = h) 

y/a 


The three pressure contributions p u , p^, and p v are defined by 

the formulas of the preceding section at various time intervals indicated 
by the digits 0, 1, 2, 3* • • • 10. ^ The time intervals were chosen as 
follows: The time interval denoted by 3 represents the end of the first 

phase of the undisturbed plate (compare section 4). At rj = 1 (period: 4) 
the influence of the adjacent plates begins and ends at t) = 0.244 
(period: 9)* At T) = 0.707 (period: 6 ), the Mach line emerging from the 
leading edge of the plate strikes the adjacent plate. 

Figure 54 represents the position of the gust front and the area 
disturbed by it at the different time intervals . 

Figure 55 illustrates the pressure of the undisturbed plate p . 

Figure 56 illustrates the pressure contribution p^. 

The pressure contribution p v is computed graphically, the velocity 
distributions Vy/vQ required for it are obtained by equation (122) for 
the time intervals 5 , 6 , 7 , 8 and reproduced in figure 57 - 


^%he corresponding curves in figures 55 > 5 ^, 57 > 60 , 6 l so:e 

denoted by the same digits. 
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Since the integrand 


Vy/VQ 



becomes infinite at the two limits and £2^ which are, as 

known, the solution of the function below the root, the graphical solu- 
tion is continued to - e) and ^£3^^ “ € 'J ^ ere e > e * axe small 

real values in comparison to 

Figure 58 represents several of the functions f for different time 
intervals . 

The integration over e and e ! is made analytically, by putting 
Vy / vq = .constant mean value. 

The relation for t^/t at rj < O.707, that is, when the additional 
source distribution is bound by a Mach line, is 

t]_/t = r| cos p. + ^ sin ^ 

If the source distribution is limited by the circle of radius R = at, 
a similar relation 

t-^/t = Tj COS P + Cl sin 3 


is applicable (fig. 59)- 

Figure 60 represents the pressure contribution p y , figure 6 1 the 
resultant pressure distributions. 

The lift of the plate (fig. 62 ) is obtained by graphical integration 
over the resultant pressures. At the appearance of the adjacent plate 
effect the lift decreases with the time interval; A 1 represents the 
steady lift of the undisturbed plate. 
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Figure 65 shows the moment distribution M plotted against plate 
center; M s t represents its steady value. Here the moment increases with 

the time because of the built-up negative pressure from t = h/a. 
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CHAPTER VII. EFFICIENCY OF A SUPERSONIC PROPELLER 
1. Introduction 


The cascade efficiency defined from thrust and tangential force is 
suitable also for the propeller. But in the preceding arguments the flow 
was assumed parallel and the blades as infinitely thin plates, which now 
must be modified. The friction at the plates must be allowed for and 
the infinitely thin plates replaced by profiles of finite thickness. 

Then the results are used to calculate the efficiency of a real propeller 
in order to obtain an approximate picture of* the efficiency to be expected. 


2. Effect of Friction on Cascade Efficiency 

When the friction at the plate surfaces is taken into account, the 
resultant force K without friction defined in chapters III and IV, is 
supplemented by an additional resistance F, so that K* is the total 
force acting on the plate (fig. 64 ). 

The frictional force is parallel to the plate. But since its com- 
ponent normal to the airstream direction is small at the angles of attack 
in question, the total frictional force can be assumed to be in the flow 
direction. 

The drag coefficient is expressed by 


c 


1 

w 


4 * 2 + gc f 

\j M 2 - 1 


( 137 ) 


Cf = is coefficient of friction of one side of the plate and q-^ 

2qi 

is dynamic pressure of inflow. 

The lift coefficient remains 



(138) 


as for parallel flow. 
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According to the definition introduced in chapter III, the cascade 
efficiency is 


T) = 


S 1 Y 
T'u 


( 139 ) 


where S' and T' are the thrust and tangential force corresponding 
to the new force K' . 

In terms of angles |3, y and a (fig- 64) the efficiency is 

H = tan a tan(p - \|f) 


where 


a = 


tan 


§; 

T' 


( 140 ) 


Introducing the drag/lift ratio 


, = W_t _F = 

A c a 



(l4l) 


the efficiency t] becomes 


1 - € t tan(p - i|/) 
1 + e 1 /tan(p - \|f) 


( 142 ) 


Hence it is apparent that, contrary to the earlier results, the 
efficiency is now dependent on the Mach number. 
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With the assumption of a turbulent boundary layer 


c f 


0.455 


iogio 



(143) 


according to Schlichting, where Re is the Reynolds number based upon 
chord l and relative velocity w. 

The values plotted in table 9 and figure 65 as functions of the 
angle of stagger were calculated with the Mach numbers M = 1.40 and 
M = 2.50, then at an angle \|r = 3° and the optimum anglesl3 i|r = 2.65° 
(M = 1.40) and \|r = 4.11° (M = 2.5O). The Reynolds number assumed at 
Re = 106 corresponds to Cf = 0.0045- 

The effect of friction is illustrated in figure 65, along with the 
efficiency curve for i|r = 3° with friction discounted. 


3. Effect of Thickness 

To assure minimum wave resistance the contour of a supersonic profile 
must consist of straight lines and its maximum thickness lie in the center, 
hence a double-wedge profile is recommended (fig. 66). 


By linearized airfoil theory (ref. 1) the thickness causes a drag 
which increases quadrat ically with the thickness ratio d/l, and which 
can be directly superposed on the lift coefficient and the frictional 
drag of the plate. 


Hence the drag coefficient of a profile of finite thickness ratio 
with friction is 


'W 


4 


M 2 - 1 


i 2 + ( - 1 + c f 


\Jm2 - 


(144) 


But the lift coefficient remains unchanged 


4\|r 



13 

As stated in the introduction, 
the plate, obviously d/l = 0. 


*opt = \ UA) 2 + c f 


\I M 2 - 


for 
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and the drag/lift ratio to be inserted in equation (l 4 l) in place of e' 
is 


!! 

G 


C W 

c a 


= \|f + ± 
♦ 


ft + 



(145) 


Table 10 shows the efficiencies of two 
profiles and the relative maximum thickness 


cascades of double-wedge 
ratio 


d/l = 0.05 and 0.10 - with friction 
at M = 1.40 

♦ = 3° 

Cj = 0.00^5 

These values are also shown in figure 67 together with those for 
d/l = 0 (the plate) for comparison. 

The angle of stagger 3 - with friction and finite thickness - 
for maximum efficiency at fixed angle of attack and fixed Mach number is 
found by simple differentiation at 



(146) 


4. Appraisal of the Efficiency of a Supersonic Propeller 

On a supersonic propeller the blades are struck at a relative speed 
which at every point of the blade is greater than the sonic velocity. 

Two types of propellers are differentiated. The one moves forward at 
supersonic speed, so that supersonic speed occurs at every rpm and every (3 . 
On the other the supersonic speed is reached without it having to move for- 
ward with supersonic speed. The efficiency of the first type propeller 
is calculated. 
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The propeller has a forward speed v of about 406 m/sec (M = 1.20 - 
sonic velocity a = 358 m/sec); it has four blades of 2m outside diameter 
and 1m inside diameter and 0.5 bub ratio. The cross section of the blade 
is a double-wedge profile, with maximum thickness ratio of d/Z = 0.07 
at the hub, and tapering to d/Z = 0.05 at the tip. 

The maximum efficiency of a profile is reached with (3 = 45°, 
according to chapter III. At the blade tip where the thrust is highest, 
this condition gives a tip speed of 400 m/sec; that corresponds to 
3, 820 rpm. For reasons of strength the blade chord tapers from = 40 cm 

at the hub to Lg = 30 cm at the tip. 

In each coaxial cylindrical section - with respect to the propeller 
axis - the angle between the relative flow direction and the profile axis - 
the angle of attack - was assumed at ^ 0 pt (compare introduction). To 

satisfy this condition, the angle of stagger in each section, that is, 
the angle between profile axis and direction of peripheral speed U must 
be varied. The relation 


tan(p - \|r) = v/u 


must be satisfied. 

With reference to a system of coordinates fixed in space, each point 
of the propeller moves on a helical line. Disturbances issue from each 
point which at the assumed pressure conditions and angles of attack can 
be regarded only as sound disturbances. The zone disturbed by each blade 
is then limited by the enveloping curve of all spheres whose centers lie 
on the various helical lines and whose radii at the same time are equal 
to sonic velocity x time. Figure 68 represents the disturbed zone of 
an edge 0A, which, for example, moves at a forward speed of 1.2 X a 
and whose maximum tip speed equals the sonic velocity; 0 ! A* represents 
the position of the same edge after a time interval At, which corresponds 
to a fourth of a revolution. 

Considering that the blades are twisted, that the disturbances of 
different sections can influence one another and be reflected on propeller 
hub, it is readily apparent that an exact calculation of the forces on 
each blade represents a difficult problem. When each blade is outside 
the zone of disturbance of the other blade, the blade can be examined 
separately. Assuming homogeneous flow and coaxial cylindrical areas, 
that is, radial equilibrium, the blade forces can be determined from a 
two-dimensional consideration of the developed blade (fig. 69 ), by com- 
puting the lift and drag and from it the thrust and tangential force in 
each cross section by linearized theory. 
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At the velocities selected the boundary effect is confined to a 
moderately large zone compared to blade area, so that its effect within 
the framework of the intended appraisal on the total forces can be 
disregarded. 

The thrust of the whole blade is then 


blade 


r ti p ds 

1 — dr 

hub 


(147) 


The integration is made by graphical method (fig. 71(a)) with 


dS 

dr 


4 \|r pw 

U - 1 


2 cos|_(p - \Jr) + 7] 


cos 7 


(148) 


computed for five sections (fig. 69). 

The corresponding Reynolds number for all sections was assumed at 
Cf = 0.004 (turbulent boundary layer). 

The torque D of a blade is defined the same way as the thrust 
by integration (fig. 7l(l))' The following relation applies: 


4^3 lade ~ 



— r dr = 
dr 



sin p - i|r) + 7] 

— r dr 

cos 7 


(149) 

The characteristics for the five sections are correlated in table 2. 
The integration gives 

s blade = 

°blade = 600 kg/m 
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hence for the propeller 


S = 1700 kg 


D = 2^00 kgm 


The efficiency of the propeller is given by 


Sv 


( 150 ) 


n = 


2jthD 


with the values inserted gives 


r| = 71.8 percent 


A quick and close estimate of the efficiency is obtainable directly 


1. Lift and drag coefficients for the flat plate at various Mach 
numbers, ranging from 1.20 to 10 , and for different angles of incidence 
are calculated, account being taken of the exact flow over both sides of 
the plate. These values are tabulated and also given in the form of. 
charts. The same coefficients are also calculated under the assumptions 
of linearized flow over the plate, according to the Ackeret theory. A 
comparison of both methods shows reasonable agreement between the lin- 
earized theory and the exact method within the usual range of angles of 
incidence (max 10°) and for the usual Mach numbers. Special formulas 
for calculating the lift and drag coefficients for very high Mach num- 
bers are derived. 

2. An analytical solution of the problem of the interaction between 
shock waves and expansion waves has been established. 

3 . A method for calculating the lift and drag coefficients for a 
cascade of flat plates is described and applied to an example, with the 
aid of the formulas derived in the foregoing item. A definition for the 


from the calculation of 



and [ — ) , where these values are appli- 

l dr L, 


cable to the whole blade. 


SUMMARY 
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efficiency of the cascade - without friction - is introduced and the 
efficiency is evaluated for two Mach numbers and different angles of 
blading. 

4. A linearized theory for supersonic flow through a cascade of flat 
plates is established and applied to the example already treated. Com- 
parison of the lift coefficients shows reasonable agreement. 

5 . For demonstration purposes, schlieren photographs were made showing 
the flow between two flat surfaces. They serve to confirm the established 
linearized theory for small angles of incidence and show clearly the inter- 
action between shock and expansion waves. 

6 . Under the assumption that the flow through the cascade of flat 
plates undergoes a small sudden change of direction, that is, a small 
change in the angle of incidence, the nonstationary flow in the cascade 
is discussed to show the kind of forces which act on the plates during 
the transition period. An example has been calculated in detail. 

7 . The definition of the efficiency mentioned in 3.> is especially 
suitable for application to a supersonic propeller. The effect of fric- 
tion and blade thickness on that efficiency is shown. A rough estimation 
of the efficiency of a supersonic propeller is then made. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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TABLE 1 


Mi 

c? 0 

Os 

?s° 

Pl/PO 

<a/po 

<i/pi 

1.00 

0 

90 

0.52830 

0.36981 

0.700 

1.10 

1.40 

73-68 

.46835 

.39704 

.847 

1.20 

3.70 

68.08 

. 41238 

.41567 

1.008 

1.30 

6.32 

65.12 

. 36092 

.42689 

I.I83 

1.40 

9-03 

63.33 

. 31424 

. 43114 

1.372 

1.50 

11.67 

62.25 

.27240 

.43050 

1-575 

1.60 

14.24 

61.65 

.23527 

.42182 

1.792 

1.70 

16.63 

61.37 

. 20259 

.40995 

2.023 

1.80 

18.84 

61.28 

.17404 

.39472 

2.268 

1-90 

20.87 

61.35 

.14924 

•37714 

2.527 

2.00 

22.71 

6 i .48 

.12780 

•35750 

2.800 

2.20 

25-90 

61.90 

.09352 

.31684 

3.388 

2.50 

29.67 

62.40 

.05853 

.25610 

4.375 

3.00 

34.01 

63.77 

. 02722 

.17143 

6.300 

4.00 

38.75 

65.25 

.00658 

•07375 

11.20 

5-00 

4 l.ll 

66.20 

.00189 

.03306 

17.50 

6.00 

42.44 

66.75 

.00063 

.01588 

25.20 

8.00 

43-79 

67.00 

.00010 

.00448 

44.80 

10.00 

44.43 

67.12 

.00002 

.00165 

70.00 

00 

45.58 

67.70 

0 

0 

00 
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TABLE 2 

LIFT AND DRAG COEFFICIENTS 


Mi 

*° 

p 2 /p 0 

7 0 

P 2 '/Pl 

p 2 7 p 0 

c a 

c v 

e 

1.20 

1 

0.39145 

58.75 

1.056 

0.43527 

0.1054 

0.0018 

0.0175 


2 

. 37210 

61.10 

1.120 

. 46187 

.2158 

•0075 

.0349 


3 

.35403 

64.37 

1.199 

.49444 

•3373 

.0177 

.0524 


3-7 

.35952 

68.08 

1.277 

.52681 

.4011 

.0259 

.06L1 

1.4o 

1 

. 29910 

46.87 

1.051 

.33027 

•0723 

.0013 

•0175 


2 

.28480 

48.19 

1.104 

. 34692 

.1440 

.0050 

.0349 


3 

. 27114 

49-57 

1.159 

. 36420 

.2156 

.0113 

.0524 


k 

.25824 

51.15 

1.219 

.38306 

.2888 

.0202 

.0699 


6 

.25376 

54.62 

1-353 

. 42517 

.4415 

.0464 

.1052 


8 

. 21130 

59.36 

1.527 

. .47984 

.6168 

.0867 

.11406 


9.03 

.20050 

63.17 

1.655 

.52007 

•7321 

.1159 

.1584 

I.60 

1 

.22374 

39.67 

1.060 

•24939 

.0608 

.0011 

•0175 


2 

.21269 

40.73 

1.104 

.25978 

.1116 

.0039 

•0349 


3 

.20207 

41.82 

1.161 

.27300 

.1679 

.0088 

.0524 


k 

.19191 

42.93 

1.219 

.28679 

.2244 

.0157 

.0699 


6 

.17280 

45.36 

1-345 

.31637 

.3385 

.0356 

.1051 


8 

.15525 

48.04 

1.484 

.34938 

•4557 

.0614-1 

.l4o6 


10 

.13914 

51-14 

1.644 

.38685 

.5783 

.1019 

.1762 


12 

. 12438 

54.89 

1.832 

. 43101 

.7110 

.1511 

.2125 


14.24 

. 11022 

61.65 

2.143 

. 50418 

.9052 

.2532 

.2532 

I.80 

1 

.16503 

34.64 

1.054 

.18352 

.0468 

.0008 

•0175 


2 

. 15640 

35-53 

1.110 

. 19318 

•0931 

.0033 

.0349 


3 

. 14813 

36.48 

1.170 

. 20363 

.1404 

.0074 

.0524 


4 

.14019 

37-44 

1.230 

. 21407 

.1867 

.0131 

.0699 


6 

.12534 

39-49 

1.362 

.23704 

.2814 

.0296 

.1051 


8 

.11175 

41.69 

I.505 

.26193 

.3768 

.0530 

.l4o6 


10 

•09935 

44.06 

I.661 

.28908 

• 4734 

.0834 

•1763 


12 

.08806 

46.70 

1.835 

.31936 

•5732 

.1218 

.2126 


15 

.07303 

51-35 

2.139 

•37227 

•7323 

.1962 

.2679 


18 

.06011 

58.00 

2.551 

. 44398 

• 9249 

.3005 

•3249 


18.84 

.05788 

61.28 

2.740 

.47687 

1.0045 

.3427 

• 3412 

2.00 

1 

.12076 

30.82 

1.058 

•13521 

.0404 

.0007 

.0175 


2 

. 11401 

31.65 

1.118 

.14288 

.0807 

.0028 

.0349 


3 

•10757 

32.58 

1.181 

1.15093 

.1211 

.063 

•0523 


4 

. 10141 

33-40 

1.247 

•15937 

.1618 

.0113 

.0699 


6 

.08994 

35-24 

1-377 

. 17726 

.2429 

.0255 

.1051 


8 

.07949 

37-22 

1-539 

. 19668 

• 3246 

.0456 

• l4o6 


10 

.07005 

39-32 

1.707 

. 21815 

.408 

.0719 

•1763 


12 

.06149 

41.59 

1.889 

. 24141 

.4923 

.1046 

.2125 


15 

.05024 

45.34 

2.195 

.28052 

.6222 

.1667 

.2679 


18 

.04070 

49.78 

2.555 

.32653 

.7604 

.214.71 

•3249 


21 

.03265 

55.67 

3-014 

.38519 

.9207 

•3534 

.3838 


22.71 

.02886 

61.48 

3-460 

. 44219 

I.0659 

.4462 

.4187 
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TABLE 2 .- Concluded 


LIFT AND DRAG COEFFICIENTS 


Mx 

+° 

P2/P0 

y o 

P 2>1 

P 2 ’/P 0 

c a 

c v 

e 

2.50 

1 

0.05296 

24.35 

1.068 

0.06251 

0.0373 

0.0006 

0.0175 


2 

.05113 

25.05 

l.l 4 l 

.06678 

.0611 

.0021 

.0349 


3 

.04772 

25.82 

1.216 

.07117 

.0914 

.0048 

.0524 


4 

.04450 

26.62 

1.296 

•07585 

.1221 

.0085 

.0699 


6 

.03857 

28.27 

1.452 

.08498 

.1826 

.0189 

.1051 


8 

.03^55 

30.00 

1.658 

.09704 

.2416 

.0340 

.1406 


10 

.02859 

31.86 

I.865 

.10916 

.3098 

•0536 

.1763 


12 

.02445 

33.81 

2.091 

.12232 

•3738 

•0795 

.2126 


13 

.01917 

36.95 

2.467 

.14439 

.4723 

.1266 

.2679 


18 

.01484 

40. 40 

2.895 

.16949 

•5742 

.1865 

• 3249 


22 

' .01035 

45.62 

2-557 

.20816 

.7162 

.2893 

• 4 o 4 o 


28 

•00575 

56.35 

4.885 

• .28592 

• 9659 

• 5136 

• 5317 


29.67 

.00482 

62.65 

5-602 

.44220 

1.0960 

.6240 

• 5695 

5-00 

3.60 

.00118 

14 

1.541 

.00291 

•0523 

•0033 

.0627 


6.17 

.00084 

16 

2.051 

.00388 

.0914 

.0098 

.1076 


10 . 68 

.00043 

20 

3-247 

• 006 l 4 

.1698 

.0319 

.1879 


16.60 

.00018 

26 

5.436 

.01027 

.2926 

.0869 

.2972 


20.21 

.00009 

30 

7-129 

.01347 

.3800 

• 1393 

.3666 


26.78 

.00002 

38 

10.893 

.02059 

•5557 

.2792 

• 5024 


31.21 

.000009 

44 

13.913 

.02629 

.6805 

.4088 

.6000 


36.29 

.000003 

52 

17-953 

•03392 

.8284 

• 6 o 48 

.7301 


4 l.ll 

.000001 

66 

24. 206 

•04575 

1.0427 

.9098 

.8726 

10.00 

3-21 


8 

2.091 


.0228 

.0013 

.0560 


7.63 


12 

4.877 


.0614 

.0082 

.1340 


13.31 


18 

10.981 


.1549 

.0365 

•2356 


18.53 


24 

19.135 


.2613 

.0877 

•3351 


23.47 


30 

29.018 


• 3828 

.1662 

.4341 


28.17 


36 

40.164 


• 5079 

.2726 

• 5366 


32.59 


42 

52.080 


.6288 

.4017 

•6393 


36.64 


48 

64.263 


• 7488 

.5565 

.7437 


40.18 


54 

76.213 


.8641 

.7292 

.8441 


42.90 


60 

87.361 


• 9529 

.8862 

•9293 


44.43 


67.12 

98.713 


1.0125 

• 9925 

.9803 

CO 

8-37 


10 



.0497 

• 0073 

.1471 


16.57 


20 



• 1873 

.0558 

.2973 


24.50 


30 



•3795 

•1731 

•4557 


28.31 


35 



.4741 

.2559 

.5386 


32.05 


40 



• 5843 

.3662 

.6261 


35.55 


45 



.6812 

.4846 

.7146 


38.81 


50 



.7840 

.6121 

.8042 


4 l .6 


55 



• 8735 

.7754 

.8876 


43.9 


60 



• 9452 

• 9123 

.9623 


45.58 


67.8 



1.0001 

1.0204 

1.0200 

N 

5 





.0152 

.0013 

•0875 


10 





.0594 

.0105 

.1763 


15 





.1294 

•0347 

.2679 


20 





.2198 

.0801 

• 3640 


25 





•3238 

.1509 

.4663 


30 





•4330 

.2500 

•5774 


35 





.5388 

.3772 

.7002 


40 

1 f— 





•6331 

.5311 

.8391 


45 





.7068 

.7068 

1.0000 


50 





•7547 

.8992 

1.1918 
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TABLE 3 

LIFT AND DRAG COEFFICIENTS OF THE BOUNDARY CURVE 


M 1 

O 

CO 

P 2 /P 0 

7 s 

0 » 

P 2 '/ P l 

P 2 / P 0 

c a 

c w 

G 

1.10 

1.4 

0.43110 

73 

14 

1.130 

0.52924 

0.2471 

0.0104 

0.0248 

1.20 

3.70 

•35952 

68 

5 

1.277 

.52681 

.4011 

.0259 

.0641 

1.30 

6.32 

.26668 

65 

7 

1.457 

.52586 

.6035 

.0668 

.1107 

l. 4 o 

9.03 

. 20050 

63 

10 

1.656 

. 52009 

•7321 

.1159 

.1584 

1.50 

II.67 

. 14392 

62 

15 

1.894 

.51593 

•8339 

.1724 

.2068 

1.60 

14.24 

.11022 

61 

39 

2.143 

. 50418 

.9052 

.2297 

.2532 

1.70 

16.63 

.07918 

61 

22 

2.439 

.49412 

.9697 

.2897 

.2988 

1.80 

18.84 

.05788 

61 

17 

2.741 

.47687 

1.0045 

.3427 

• 3412 

1.90 

20 . 87 

.04047 

61 

21 

3.097 

.46220 

1.0449 

.3983 

.3812 

2.00 

22.71 

.02886 

61 

29 

3.460 

.44219 

1.0665 

.4465 

.4187 

2.20 

25.90 

.01417 

61 

54 

4.250 

.38746 

1.0883 

.5284 

.4855 

2.50 

29.67 

.00482 

62 

39 

5.602 

. 44221 

1.0960 

.6242 

.5695 

3.00 

34.01 

.00073 

63 

46 

8.385 

.22824 

1.0735 

.7243 

.6747 

4.00 

38.75 

.000014 

65 

15 

15.25 

. 10060 

1.0650 

.8551 

.8025 

5.00 

41.11 

.000001 

66 

12 

24.40 

.04612 

I.0515 

.9107 

.8652 

6.00 

42.44 

0 

66 

45 

35.25 

.02231 

1.0362 

.9432 

.9136 

8.00 

43-79 

0 

67 

0 

63.ll 

.00630 

1.0171 

.9692 

.9601 

10.00 

44.43 

0 

67 

7 

98.90 

.00233 

1.0121 

.9895 

.9803 

00 

45.58 

0 

67 

41 



1.0050 

1.0185 

1.0200 
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TABLE 4 ' 


LIFT AND DRAG COEFFICIENTS ACCORDING TO THE LINEARIZED THEORY 


Mi 


c a 

c v 

€ 

1.20 

v 1 

0.10353 

0.00018 

0.0175 


2 

.21046 

.00734 

.0349 


3 

.31600 

.01655 

.0524 


3.7 

•30957 

.02516 

.0646 

1.40 

1 

.07144 

.00125 

.0175 


2 

. 14248 

.00497 

.0349 


3 

.21392 

.01121 

.0524 


4 

.28496 

.01989 

.0698 


6 

.42744 

.04480 

.1047 


8 

.56992 

•07950 

.1396 


9.03 

.64259 

.10112 

•1574 

1.60 

1 

.05604 

.00089 

•0175 


2 

.11177 

.00390 

.0349 


3 

. 16781 

.00879 

.0524 


4 

.22353 

.01560 

.0698 


6 

•33530 

.03458 

.1047 


8 

.44707 

.06240 

.1396 


10 

.55884 

.09752 

.1745 


12 

.67060 

.14040 

.2094 


14.24 

•79582 

•19772 

.2485 

1.80 

1 

.04676 

.00082 

•0175 


2 

.09325 

.00326 

.0349 


3 

.14001 

.00734 

.0524 


4 

.18650 

.01302 

.0698 


6 

.27976 

.02932 

• 1047 . 


8 

•37301 

•05205 

.1396 


10 

. 46626 

.08135 

.1745 


12 

•55952 

.11714 

.2094 


15 

.69953 

.18314 

.2618 


18.84 

.87855 

.28884 

.3288 

2.00 

1 

.04042 

.00070 

•0175 


2 

.08060 

.00281 

.0349 


3 

.12102 

.00634 

.0524 


4 

.16120 

.01125 

.0698 


6 

.24180 

.02532 

.1047 


8 

•32240 

.04995 

.1396 


10 

.40301 

.07030 

.1745 


12 

.48361 

.10125 

.2094 


15 

.60463 

.15829 

.2618 


18 

. 72564 

.22802 

.3142 


21 

.84643 

•31040 

.3665 


22.71 

.91502 

•36259 

.3962 

2.50 

1 

.03054 

.00053 

•0175 


2 

.06091 

•00213 

.0349 


3 

.09145 

.00479 

.0524 


4 

.12181 

.00851 

.0698 


6 

.18272 

.01915 

.1047 


8 

•24363 

.03400 

.1396 


10 

.30454 

.05318 

.1745 


12 

.36544 

.07650 

.2094 


15 

.45689 

.11962 

.2618 


18 

.54834 

.17230 

• 3142 


22 

.67016 

.25742 

• 3840 


28 

.85288 

.41675 

.4887 


29.67 

.90366 

.46789 

.5178 

5-00 

4 

.05699 

.00398 

.0698 


6 

.08549 

.00896 

.1047 


12 

.17097 

.03579 

.2094 


18 

.25654 

.08061 

• 3142 


24 

•34203 

.14321 

•4189 


30 

.42752 

•2?385 

.4712 


36 

.47030 

.32232 

.5760 


41.11 

.58584 

.42033 

•7175 

10.00 

4 

.02806 

.00196 

.0698 


6 

•04209 

.00441 

.1047 


12 

.08418 

.01762 

.2094 


18 

.12631 

.03969 

• 31^2 


24 

.16840 

.07051 

.4189 


30 

.21049 

.11022 

.4712 


36 

.25258 

.15800 

.5760 


42 

.29467 

.21600 

•7330 


44.43 

.31172 

.24177 

•7755 
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TABLE 5 


Mi 


( c aG " c al) 

C aG 

( c vG ~ c wl) 

c wG 

1.40 

1 

0.00084 

0.07228 

0 

0.0013 


2 

.00151 

•14399 

0 

.0050 


3 

.00163 

.21555 

.0001 

.0113 


4 

.00385 

.28881 

.0002 

.0202 


6 

.01408 

.44152 

.0010 

.0464 


8* 

.04690 

.61682 

.0072 

.0867 


9-03 

.08950 

.73209 

.0148 

•1159 

5.00 

4 

.0005 

.0523 

.0009 

.0033 


6 

.0030 

.0914 

.0010 

.0098 


12 

.0255 

.1698 

.0065 

.0319 


18 

.0696 

.2996 

.0246 

.0869 


24 

• 1391 

.3800 

.0341 

•1393 


30 

.1846 

•5557 

• 1453 

.2792 


36 

•3972 

.6805 

.2698 

.4088 


41.11 

. 4669 

.8284 

.4895 

.6048 
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TABLE 6 


Region 

P/POi 

M 

M-° 

PO/PO! 

1 

0 • 31404 

l.4oo4 

45.5 

1 

2 

.3645 

1.293 

50.6 

.9996 

3 

• 2991 

1.435 

44.2 

1 

4 

•3455 

1.332 

48.8 

• 9997 

5 

.2848 

1.469 

42.9 

l 

6 

• 3304 

1.364 

47.1 

• 9 997 

7 

.2711 

1.503 

41.7 

l 

8 

.3139 

1.401 

45.7 

.9998 

9 

.3285 

1.368 

47.0 

• 9997 

10 

• 3139 

1.401 

45.6 

• 9997 

11 

.2985 

1.436 

44.1 

. 9997 

12 

.2999 

1.433 

44.2 

• 9997 

13 

.2850 

1.468 

42.9 

• 9997 

i4 

.2708 

1.504 

41.7 

■ 999 7 

15 

.3624 

1.154 

50.1 

.99 82 

16 

.3445 

1.334 

48.6 

■ 9991 

17 

• 3295 

1.366 

47.0 

■ 999 7 

18 

.3124 

1.404 

45.4 

■ 9991 

19 

.3276 

1-370 

46.9 

■ 9998 

20 

.3134 

1.402 

45.5 

• 9997 

21 

.2972 

1.439 

44.0 

■ 9991 

22 

.2996 

1.434 

44.8 

■ 9991 

23 

.2842 

1.483 

42.8 

■ 9991 

24 

25 

.2696 

.3610 

1.507 

41.6 

■ 9991 
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TABLE 7 


(a) Angle of deflection and shock angle of the 
shocks between the zones 


Regions 

Angle of ^.deflection 
5° 

Shock angle 
7° 

1-2 

3-00 

49-57 

3-4 

5-03 

47.91 

5-6 

2.93 

46.48 

7-8 

2.93 

45.13 

8-15 

2.93 

49.43 

11-16 

2.92 

47.78 

13-17 

2.90 

46.40 

14-18 

2.93 

45.03 

18-25 

2.93 

49.30 


(b) Intensity of expansion waves 
between the zones 


Regions 

Intens ity 
£0° 

Regions 

Intensity 

£8° 

1-3 

1.00 

13-14 

1.02 

2-4 

1.03 

15-16 

1.00 

5-5 

1.00 

16-17 

.88 

4-6 

.89 

16-19 

1.00 

4-9 

1.03 

17-18 

1.06 

5-7 

1.00 

17-20 

• 99 

6-8 

1.00 

18-21 

1.00 

8-11 

1.01 

19-20 

.89 

9-10 

• 92 

20-21 

1.05 

10-11 

1.00 

21-23 

.89 

10-12 

• 92 

22-23 

1.05 

11-13 

• 9i 

23-24 

1.05 

12-13 

1.02 
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TABLE 8 


CASCADE EFFICIENCY i\ PERCENT = f(p,+) (NO FRICTION) 


( 3 ° 

1 

3 

5 

7 

10 

10 

89.8 

69.6 

49.6 

29.72 


20 

94.7 

83-9 

73.6 

63.4 

48.4 

50 

96.0 

88.3 

80.8 

73-5 

63.1 

4-0 

96.5 

89.8 

83.4 

77-3 

68.8 

50 

96.6 

90.00 

8 3.8 

78.2 

70.4 

60 

96.0 

88.8 

82.4 

76.5 

68.8 

70 

94.8 

85.7 

78.2 

71.4 

63.1 

80 

90.6 

76.4 

65.8 

57.65 

48.4 

3° for max. t) 

45-5 

46.5 

47.5 

48.5 

50.0 

%iax> percent 

96.3 

90.2 

83.7 

78.2 

70.4 
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TABLE 9 


CASCADE EFFICIENCY WITH FRICTION ALLOWED FOR 


3° 

M = 

i.4o 

M = 

2.50 

i = 3° 

t = 2.65 

+ = 3° 

i = 4.11 


10 

55.8 

57-0 

44.1 

4l.l 

20 

74.2 

74.6 

63.9 

63.7 

30 

80.3 

80. k 

71.2 

71.8 

4o 

82.5 

82.8 

73-9 

74.8 

50 

8 2.6 

82.7 

73-5 

74.7 

60 

80.5 

80.6 

69.9 

71.8 

70 

-74-7 

74.6 

60.6 

63.8 

80 

57-8 

57-2 

33-5 

4l.4 


TABLE 10 


CASCADE EFFICIENCY WITH DOUBLE-WEDGE PROFILE AND FRICTION 

M = 1.40 = 3° 


0° 

d jl = 0.05 

d/l = 0.10 

10 

!+5-53 

29.01+ 

20 

65.28 

47.22 

30 

72.52 

74.78 

l+o 

75-11 

56.93 

50 

7^.83 

54.82 

6o 

71.50 

47.31 

TO 

62.71 

29.25 

8o 

37-18 
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TABLE 11 


Section 

1 

hub 

N 

2 

3 

center 

M 

4 

5 

tip 

S 

Radius, rm 

0.500 

0.625 

0.750 

0.875 

1.000 

Profile chord, m 

0.450 

0.425 

o.4oo 

0.375 

0.350 

Thickness ratio, d /l 

0.070 

0.065 

0.06 

0.055 

0.050 

Max. thickness, dm 

0.0315 

0.0282 

0.024 

0.0206 

0.0175 

Tip speed, um/sec 

200 

250 

300 

350 

400 

Relative velocity, wm/sec .... 

L52.7 

476.8 

504.6 

536.0 

570.0 

Relative Mach number, M 

1.339 

1.411 

1.493 

1.586 

1.687 

Angle of attack, ^ont* 4 £ • • • • 

0.0818 

0.0787 

0.0761 

0.0742 

0.0721 

opt | Q 

4.6 7 

4.52 

4.38, 

4.25 

4.13 

Lift coefficient, c a 

0.3672 

0.3164 

0.2745 

0.2407 

0.2120 

Drag/lift ratio, € 

0.1637 

0.1572 

0.1523 

0.1484 

0.1442 

Gliding angle, 7 0 

9.30 

8.95 

8.66 

8.43 

8.21 

Stagger, p° 

68.8 

62.9 

57.9 

53.5 

49.58 

dS/dr, kg/m 

612 

748 

825 

878 

904 

dD/dr, kg 

1026 

1118 

1173 

1212 

1224 
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Figure 1.- Expansion around a corner. 



Figure 2. - Oblique compression shock. 
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Figure 3.- Shock and flow deflection angles for Mg = 1. 
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Figure 4.- Flat plate in supersonic flow. 



0 1 23456789 10 M 1 

(a) Dynamic pressure/stagnation pressure 

(b) Dynamic pressure/ inflow pressure 


Figure 5. 
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Figure 7 .- Drag of flat plate. 
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Figure 8 .- Polars of flat plate. 
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Figure 9.- Flat plate in linearized flow. 



Figure 10.- Lift by linearized theory. 
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Figure 11.- Drag by linearized theory. 


Characteristics 
Shock polar 






(a) Overtaking 

(b) Crossing (in free jets) 

(c) Reflection 

(d) Cascade flow 

Figure 12 .- Interaction between expansion waves and compression 

shocks. 
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Figure 13.- Small expansion. 



Figure 14. - Small variation for one compression shock. 



Figure 15.- Compression shock overtaken by expansion wave. 


14K 


MCA 1M 1369 


101 



Figure 16.- Schematic representation of the graphical solution. 
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Figure 18.- Expansion wave crosses compression shock. 
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Figure 19.- Schematic representation of graphical solution. 



Figure 20. - Crossing of expansion waves. 
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Figure 22.- The cascade problem. 



Figure 23.- Cascade example. 
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Figure 24.- Directional changes of interference lines of the cascade 

example. Large digits denote zones; small digits denote directional o 

changes at crossing of expansion waves and compression shocks. 01 

At reflections they indicate the direction of the reflected waves; flow 
always from left to right. 
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P/P 0 

0.6 

04 

0.2 

0 

0.2 

0 0.1 0.2 0.3 0.4 0-5 0.6 0.7 0.8 0.9 1.0 

1 / L 

Pressure at upper side 

Pressure at lower side 

/////// Resultant pressure 

Resultont pressure by iinearized theory 

Figure 25.- Pressure variation along the plate (exact method). 
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Figure 26.- Definition of thrust and tangential force. 
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Figure 27.- Cascade efficiency (no friction). 



Figure 28.- Linearization of crossing. 
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Figure 29. - Linearized cascade theory. 
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Figure 30.- Numerical example of linearized cascade theory. 
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Figure 31.- Pressure distribution over plate by linearized cascade 

theory. 
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Figure 32 


The profile 
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Figure 33. - Profile in test section of tunnel. 
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Figure 34. - Schematic representation of the bend in the upper tunnel 

wall. 
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Figure 35. - Position of optical axis and shadow formation, x = 8 cm; 

y = 4.5 cm; Z = 300 cm. 


Figure 36.- Profile at starting. 




Figure 37.- Schlieren diaphragm vertical, t = 0°, 
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Figure 40. - Photograph of crossing. 



Figure 41.- Perspective distortion of figure 23. For comparison with 
schlieren photograph in figure 39. (O represents the position of 
the light source; the finer lines represent the shadow boundaries of 
the plate and interference lines.) 
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Figure 42.- Sudden change in angle of attack. 



Figure 43.- The unsteady source. 



Figure 44.- The periodically created source distribution. 
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Figure 46.- Pressure distribution at y = constant lines. 
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Figure 47.- Velocity distribution at y = constant lines. 
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Figure 48. - The flat plate in a gust. 
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Figure 49.- Pressure distribution along the plate. 
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Figure 51.- Start of process. 
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Figure 53.- Second phase of additional pressure. 
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Figure 55. - Pressure of undisturbed plate p u . 



Figure 56.- Pressure contribution p^. 
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Figure 57.- Velocity distributions v y /v 0 at times 5, 6, 7, and 8. 
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Figure 58. - Functions f for various time intervals. 
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Figure 60.- Pressure contribution p v> 
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Figure 62. - Lift of plate. 



Figure 63. - Moment distribution M referred to plate center. 
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Figure 64.- Allowance for friction. 



(a) 

(b) With friction. M = 1.40, 
y = 3°. 

(b') With friction. M = 1.40, 
* = 2.65°. 


No friction. 

(c) With friction. M = 2.50, 
t= 3°. 

(c') With friction. M = 2.50, 
♦ = 4.11°. 


Figure 65.- Cascade efficiency. 
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Figure 66.- Allowance of finite thickness ratio. 
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Figure 67.- Efficiency of cascade of double -wedge profiles. 
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Figure 68. - Area disturbed by leading edge of a blade assumed as 

flat plate. 
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Figure 69.- Blade chord, maximum thickness, stagger, and location 

of cross sections. 


/ 


150 


NACA TM 1369 



Figure 70.- Blade form, blade form without angular rotation. 
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Figure 71. 


NACA - Langley Field, Va. 


